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Slow Broad Side Motion of a Flat Plate in a Viscous Liquid 


By SuBHAS CHANDRA Gupta, Agra, India!) 


1. Introduction 


The problem of slow broad-side motion of a flat plate in viscous liquid has 
been solved in terms of a function m which turns out to be the electrostatic 
potential of the plate kept at a constant potential in vacuum. The hydro- 
dynamical problem thus reduces to a problem in electrostatics with known 


boundary conditions. The function satisfies LAPLACE’s equation and the- 


boundary conditions can be expressed in terms of gy. The problem corresponding 
to a disc of any arbitrary shape can be easily worked out Pt ors up an 
experiment based on this analogy. 

The method has been applied to the case of a circular disc by introducing 


transforms and solving the problem as a boundary value problem in electro- . 
statics. This problem has been also worked out by Ray?) by using an ad-hoc . 
method, but the present method gives a systematic and logical approach to» 


problems of this kind. It has a distinct advantage over Ray’s method in which 
a correct form of the solution has to be guessed at the outset and later integral 
solutions have to be constructed which give desired discontinuity at the edges 
of the plate and this quite often is not easy. 


2. The Equations of Motion, Their Solution and the Electrical Analogy 


The hydrodynamical equations of slow viscous motion after neglecting the 


quadratic terms of inertia are: 


Op 2% . 
7 i eS UViw; (1) 
ee 2 
Oy BV", (2) 
Ops? = 

‘s =pV?w, (3) 

where 
. 0 02 02 


1) Department of Mathematics, Agra College. 
2) M. Ray, Phil. Mag. 21 (7), 553-558 (1936). 
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and the equation of continuity is 


Ou Ov Ow 4 
Ox et Oy ' 02 (4) 
Let ; 
, OP =, 9” Wee Lia die 
UA me) en ee ie Cae (5) 


Equation (4) now reduces to 


8p: OP 08 6 
oe Shae or ee (6) 


From equations (1), (2) and (3) we find that 
Pet op 
p= byt 2u 2, (7) 


p, being a constant. The plate coincides with the plane z = 0 and is bounded 
by the curve z = 0, f(x, y) = 0 and is moving in the direction of the z-axis with 
a small velocity V. The boundary conditions are 


=== 0. lV 
on the plate, i.e., within 
a= Os F(a) = On 
and 
T= i= ipa (0). 


at infinity. Thus the conditions expressed in terms of @ are 


a wnad 0 
on the plate, i.e., within ' 
2= 05, flay) = O08 


and must tend to zero at infinity in a manner such that 


may also tend to zero. Also we see from (5) that on the plane z = 0, u=v=0, 
i. e., there is no flow in the plane z = 0. The flow of the fluid is normal to the 


plane z= 0, and so the pressure throughout this plane must be constant. 
Therefore from (7) 


Op ; 
5, Oe, Ole 105, (8) 
outside the plate. Thus the problem may be briefly restated as follows: 


V?o=0, rite eee) 
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g=—V, within z= 07. f%, y)=0, (10) 
Op ‘ 
> ie QO} outsides g= Ue fix, y) =O, (11) 
and 
» 9P Op P 
Np praee Fries a ar and 


must all tend to zero at infinity. 

This problem is same as that of determining the electrostatic potential of 
a conducting electrified disc in the plane z = 0, bounded by the curve f(x, vy) =0, 
kept at a constant potential —V in vacuum. 


3. The Case of a Circular Disc 


As an example of the above general treatment we take the case of a circular 
disc. Without loss of generality we take the radius of the disc to be unity and 
z-axis centrally perpendicular to it. Transforming the equation (9) in cylindrical 
coordinates (7, 6, z) we get 

0? 1 0 0 
gee ae (12) 


Y OY 


Here we have omitted the term 07/06? due to symmetry about the z-axis. 
The boundary conditions are 


ga’ (C27 <1), | (13) 
op POR 2 = 0). 
Introducing Hankel transform defined by the relation 
a) =[er hierar, (15) 
0 
we find that 
oO 2 : z 
[02x Ibn) a = [ae — P| 6). 
0 
Equation (12) becomes 
2m — 
See pe). (16) 


Sr 
Since the boundary conditions (13) and (14) are of a mixed type, we invert p 


for g before satisfying them. aa 
In view of the symmetry of the problem it will be sufficient to consider the 
region z > 0. Since the potential must vanish at infinity, the appropriate 
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solution of (16) is 
p(p) =A(p)e"?* (2 > 0), (17) 


A(p) being an unknown function of #. 
Inversion for gives 


p= | Al) b Ilr P) ab. (18) 


If we insert (18) in the boundary conditions (13) and (14) we get the following 
dual integral equations: 


[PA®) hird) 4p=-V OSr<1), (19) 
[87 A(®) Jolrd) dp =0 (7 > 1). (20) 


0 


TitcHMARSH®) and BusBriIDGE’) have considered dual integral equations of 
this type and their solution is 


A(p) = —= x $e sinp . (21) 
Thus the required solution is 
Bi r =e 1 
p=—— fe Jord) “SF ap. (22) 


In this particular case Legendre transform can also be used employing 
oblate spheroidal coordinates. The result can be expressed in a neat form. 
Suppose (§, y) are oblate spheroidal coordinates related to the cylindrical 
coordinates (7, z) by 

z=én, 


and 


=(1— e991 + gH 


Transforming LAPLACE’s equation into oblate spheroidal coordinates and 
employing Legendre transform of » defined by 


_ E, C. Tircumarsu, Theory of Fourier Integrals (Oxford 1937), p. 335. 
) I. W. Buspripce, Proc. Lond. math. Soc. 44 (2), 115-129 (1938). 
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we find that®) oy 
» = ——— cot7'n, 
¢ = n (23) 


This solution can be shown to agree exactly with the more complicated form 
given by (22). 

The resistance of the fluid to the motion of the disc calculated from either 
(22) or (23) is 16 Vu. If the radius of the disc be a the resistance can be shown 
to be 16 Vayu, which is identical with that obtained by Ray?). 

I thank Dr. S. D. Nica for suggesting this problem to me. 


Zusammenfassung 


Das Problem der langsamen Querbewegung einer ebenen Platte in einer zihen 
Fliissigkeit ist zuriickgefiihrt worden auf ein Problem der Elektrostatik. Die 
Geschwindigkeitskomponenten und der Druck kénnen berechnet werden aus 
einer Funktion g, die als die gleiche erscheint wie das elektrostatische Potential 
einer Platte unter konstantem Potential im Vakuum. Die Strémungsverhiltnisse 
an einer Platte vom allgemeinen Umriss kénnen mit Hilfe der Experimente auf 
Grund dieser Analogie bestimmt werden. Der Fall einer kreisrunden Platte wurde 
analytisch behandelt mit Hilfe der Abbildungen. 


(Received: July 3, 1956.) 


Charge and Discharge of a Non-Linear Condenser 
Through a Linear Non-Dissipative Inductance 


By AsoxenpDU Mozumpb_Er, Kharagpur, India?) 


Introduction 


The capacity of a condenser with semi-conducting di-electric depends upon 
the voltage across it. This non-linearity arises from the processes of diffusion 
and recombination which occur within the di-electric, and is important at 
voltages below the breakdown value. At higher voltages, normal metallic con- 
duction takes place and the non-linearity becomes inappreciable. Considering 
the building up of space charge at a blocking electrode by free charge carriers 
MAcDONALD [1]?) has deduced the formulae 


. sinha V 3) 
aV 


5) C. J. Tranter, Integral Transforms in Mathematical Physics (Methuen, 1951), p. 100. 

1) Indian Institute of Technology, Department of Physics and Meteorology. 

2) Numbers in brackets refer to References, page 279. 
3) Note that in [1] the differential capacity (dQ/dV) is deduced, whereas here the corresponding 
integral capacity (Q/V) is used since it is the latter capacity which appears in KrRCcHHOFF’s circuit 
equation (cf. [6]). 


C=C, (1a) 
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where C is the capacity of the condenser with voltage V across it, Co the 
capacity of the condenser with vanishingly small voltage across it, x a universal 
constant (e/2 KT or e/4 KT according as one or both electrodes block), ¢ the 
electronic charge, K BorTzMANN’s constant and T the absolute temperature. 

Building up of space charge at a blocking electrode occurs in the case of 
suitably biased rectifying electrodes [2], photoconductors [3], electrolytes [4], 
transistors [5], liquid di-electrics [4] and the like, which, therefore, show strong 
nonlinear capacities up to a certain voltage limited either by di-electric break- 
down or by field emission inside the di-electric or by both. Although, in most 
physical systems, the range of voltage over which the effect is marked is quite 
small, an analysis of circuits containing such non-linear capacities is important 
in certain switching circuits and di-electric amplifiers and may be of use in the 
design of waveform generators. 

MACDONALD and BRrACHMAN [6] have studied the charge and discharge of 
such a condenser through a linear resistance in a dc circuit. They have also 
considered an ‘exponential’ condenser for which 


C= Cyexpiav| (1b) 


with the same significance of the symbols as in equation (la), and have found 
that the two solutions have the same general characteristics. An exponential 
condenser can be physically realized with the help of a very high gain amplifier, 
a condenser feedback device and a circuit whose output is proportional to the 
exponential function of the input [6]. 

In the following, we shall discuss the oscillatory charge and discharge of 
an exponential (1b) and a sinh-type-condenser (1a) through a linear lossless 
inductance in a dc circuit and obtain the voltage and current waveforms. 


The Charging of the Condenser 


Ae: 


iE 


\| po Cu: 


DPeaenst 
Figure 1 
The charging circuit. 


The differential equation of the circuit shown in Figure 1, for the process 
of charging, is 


dQ 
Lye tV=Vp, (2) 
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where Vj is the dc voltage of the source, V the voltage across the condenser at 
any instant of time ¢, L the inductance in the circuit and Q the charge on the 
condenser plates at the same instant of time ¢. 

For the exponential condenser (1b) we have 


0=C, Ver”! (3) 
and for the sinh-type condenser (1a) 


iE (<2) sinha V. (4) 


Case 1: The Exponential Condenser 


With the initial conditions at ¢= 0, V = 0 and current 7 = 0 (and hence 
dV |dt = 0), V is always positive and hence V can be written for |V| in equation 
(3). Differenting (3) twice with respect to time and using equation (2), we get 


@V . a(2+aV) (ey ede NM (5) 
dt? L+aV \ di EGF (vee Vy 
This equation can be linearized by the substitution 
dV\2 
z= (ie) 
to obtain ; 
og go eet T= 2(V,—V) (6) 
dV 1+ al ECre (i av) 
the complete solution of which is 
cl ez 2 e7uV ™ 1 S| 4 
= - / = a 
o> aie t LCi + av [V% ig eer ah (7a) 


where / is the constant of integration. In view of the initial condition at p=, 
dV |dt = 0 (i. e., Z = 0), A becomes 2/(«? L Cy) and we have, from (7a), 


1 ere 2 ae 1 
7 ae 


+e” @2VV,— 02? V2-1ltaV)? 2, (7b) 


We shall see later that V is a periodic function of time. Starting from ¢ = 0, 
therefore, the right-hand side of (7b) is to be taken with the positive sign until 
the first maximum is reached and then with the negative sign until the next 
minimum and so on. 

Integrating (7b), we obtain 


Vv 
van fet (aryl Vy aH evi Tre Vr eV. 


0 


2 


f= 
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This can be normalized to 


t 


tafeaeaireen-#-14apt, 6 
0 
where 


r= Eo, = 0 V GEOR ole (9) 


Numerical integration of (8) gives the waveform of the voltage across the 
condenser for a particular value of %). 
From (3) and (9), the normalized expression for current 1 = dQ/dt becomes 


j= bij +e wag = Dep (10) 
where 
: al C 
19 — Ch, L v . 


The positive and negative signs in (10) are to be alternately used in the same 
way as in equation (7b). 
Case II: The sinh-T ype Condenser 
In this case we have from equations (2) and (4) 


ay 
at? 


Ves 
iC 


+ tanhe V (=) 2 V secha V , (11) 


which can be put into a linear form by the substitution Z = (dV /dt)?, to give 


dZ 2 
av + 2o(tanhaV) Z = Tire. (V, — V) secha V , 
with the complete solution 
Z=Asech?a V + ee hey [cosha V + « (Vj — V) sinha V]. (12) 


oe 
With the initial conditions at t= 0, V = 0 and dV /dt = 0, i. e., Z = 0, we get 
2 
pohek rics 
and 


dV 2 i 
a= | == =, secha V [(cosha V — » +a(V)— V) sinha V}?. (13) 


The positive and negative signs are to be treated in the same way as in (7b). 
Equation (13) gives, on normalization, the charging waveform, starting with 
t = 0, for the voltage across the condenser. 
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x 
t 
= = [ cosh x [(cosh x — 1) + (% 9 — x) sinh x)71? dx , (14) 
0 


where t, x and x) have the same significances as in (9). 


The normalized charging current through the sinh-type condenser is given by 
(4) and (13) as__ 
t = +%[(coshx — 1) + (x) — x) sinh x]? , (15) 


with the same value of 7) as in equation (10). 
The signs in the expression (15) for current 7 are to be used in the same 
way as in (10). 


The Discharging of the Condenser 


, C 
Eee art 
Figure 2 


The discharging circuit. 


KIRCHHOFF’S equation of the circuit of Figure 2 for the case of discharge 


gives d20 


with usual significances of the symbols. 


Case 1: The Exponential Condenser 
During discharge, V is capable of assuming both signs. Hence using (3), 
(16) can be written as 


(av d? |V| 
Fre eee 


dV a\V| ad |V|\2) = 
+ 2a Far (a) + V =0, (16a) 


Ce eX lV) 


where | V | is the absolute value of V. This equation remains unchanged if 
V is replaced by —V and it will be sufficient for our purpose if we deal with 
the magnitude of V only. The solution of (16a) can then be obtained by the 
method previously outlined. With the initial conditions at t= 0, V = Vy and 
dV |dt = 0 (as the current is zero initially), we get 


aVe i: ste Tt 2 ee 
di ne ae ee (doe Oe) (16b) 


x [exVe («2 Ve + t= Vo) = et (a2 V2 — lo V)P? i) 
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which can be normalized to 


dx i) eae 
dt Ses Sie 


x 72, > 7 a 1/2 
Tica [e* (42 +1 — 2%) — & (x? +1 =H)", (17) 
where x, x) and t have the same significance as in (9) and the positive and 
negative signs are to be used as before. Differentiating equation (17), we get 


d?x Me Ze, -% £ (a! tye 
Gt ae le (46-1 = hy) Yee 
d Det? LN Al ade x) 
x ales) - a (1 4 4/2 
Hencerat 7 0) i.6., Y= Xo, 
Gey 1 alsa (2 + x») 
dt® |t-0 pee ORE PRA) 
and, therefore 
ax 
di* |t=0 


is negative. Thus at ¢ = 0, the voltage has a maximum and it falls with increas- 
ing time until a minimum is reached. Taking, therefore, the negative sign for 
dx/dt and using the initial condition, at ¢ = 0, x = x», we get 


t 
oe =e) a) (18) 
where 
I(x) = : e* (1 + x) [e% (x8 + 1 — 1%) — ec” (w® + 1— x)? dx (18a) 
0 
and 


I(x) = ii e* (1 + x) [e% (x3 4+ 1 — xq) —e* (x2 +1— x) 2 dx. (18b) 
0 


The above expressions are valid from time ¢ = 0 until the next minimum for x 
is attained. However, as we shall see later, x is periodic in time; so the above 
treatment remains valid for any time between a maximum and the next 
consecutive minimum for x. Equation (18), on numerical integration, gives the | 
discharge shape of the condenser voltage in the appropriate range. 

The normalized discharge current through the exponential condenser is, 


from (3) and (17) 
t= +45 [e (x8 + 1 — x) — &” (x7 + 1— x)" (19) 
with the same value of 7) as before and x considered in magnitude only. 


After ¢ = 0 and before reaching the next minimum for x, dx/dt is negative. 
Therefore, the current is also negative in that range and it becomes positive 
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only after the voltage crosses its minimum. This is because the direction of 
current through the condenser during discharge is opposite to that during 
charge (Figure 3) at ¢ = 0. 


jA 
aw / 
at discharge at charge 
Figure 3 
Diagram showing the directions of current at ¢ = 0 for charge and discharge. 


Case I1:, The sinh-T ype Condenser 


The normalized expression for dx/dt in the case of the sinh-type condenser 
can be obtained from (4) and (16). With the initial conditions, at ¢= 0, V = Vy 
(i. €., 4 = %») and dV /dt = 0 (i. e., dx/dt = 0), this becomes 

dx 


ap = ET sech~x [(cosh x — cosh x) — (x sinhx — x»sinhx,))"” (20) 


where x, *) and t have the usual significance. Differentiating (20), we get 


7% 


qe = +t sech «x tanh x [(cosh x — cosh %q) — (x sinh x — x, sinh x,)|"” 


(=) t? x sech x 
SN de 2 


Hence at ¢ = 0, 


ahd = pag x, sech x 
at? |\5=0 ew o 
and, therefore 
a2x% 
dt? |\t=0 


is negative. 

Thus at ¢ = 0, the voltage has a maximum and it falls with increasing time 
till a minimum is reached. Taking, therefore, the negative value of dx/dt and 
using the initial conditions, we get 


<= I(x) — I(x), (21) 


where 


I (%) = [ cosh x [(cosh x — cosh %) — (x sinh x — % sinh meV ax 


and 


I(x) = [ cosh x [(cosh x — cosh %») — (x sinhx — % sinh ae aX. 
0 
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The above expressions are valid from time ¢ = 0 till the next minimum for x 
is reached. Here again it will be shown that x is periodic in ¢. So that the 
above treatment remains valid for any time between a maximum and the next 
minimum for x. Equation (21), when numerically integrated, gives the shape 
of the condenser voltage for discharge in the appropriate range. 

The normalized discharge current through the sinh-type condenser is given 
by (4) and (20) as follows: 


i = +1, [(cosh x — cosh %) — (« sinh x — %» sinh %9)|'” (22) 


where 7, has the same value as before. 


Discussion on the Charging of the Condenser 


A close examination of the integrands in equations (8) and (14) shows that 
they are singular at x= 0 and for another value of « = x, defined for the 
exponential condenser by, 


1+ 6 (%, %) — #2? -—1+%,) =0 (23) 
and for the sinh-type condenser by, 
(cosh x, — 1) + (xp — 4) sinhx, =0. | (24) 


It is, therefore, necessary to examine the convergence of the integrals in 
equations (8) and (14) at x = 0 and at x = x,. Let ®(x) denote the integrand 
in any of the two equations (8) and (14). Then near the singularity at x = 0, 
the integrand in both cases can be written as 


Oxy) ax™ (a+b 4c ye --) 


where a, b,c... are constants and x is small but positive. Hence f D(x) ax 
exists at the lower limit x = 0. In a similar fashion, it can be shown that near 
the singularity at « = x,, we can write in both cases 


O(i2)=z2? Ct Bztyz24+--:) 


where again «, B,y,... are constants and z= x, — x is small but positive. 
Hence the integral i (x) dx also exists at the upper limit, viz., at x = x,. The 
importance of the existence of J(x,) lies in the fact, that if I(,) were infinite, the 
charging would be aperiodic. The above treatment shows that such is not the case. 

It is easily seen from equations (23) and (24) that for « greater than 4%, 
@(x) is imaginary. So, physically, there can be no normalized voltage greater 
than x,. Hence x, is the maximum voltage that can develop across the condenser 
in the charging process. This is further substantiated by the fact that at x = Xi; 


ax 


rama 
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in both cases [equations (7b) and (13)]; further 


ax ye eae: 
a +x) > 4) owe 
for the exponential condenser and 
2 
= = hy (sech x) (x9 — 4) (25) 


for the sinh-type condenser, k,, ky being constants. We shall eventually see 
that x, is always greater than x9. Thus at « = x,, d?x/dé? is negative in both 
cases, Showing a maximum at * = *,. 

The value of x, depends on the line voltage [equations (23) and (24)]. This 
has been calculated by the Newton-Raphson method for some representative 
values of the line voltage, and the results are tabulated in Table 1. 


Table 1 
wy 

Xo : ; 
Exponential sinh-Type 
1 1-477 1-688 
4 4°793 © 4-987 
10 10-909 11-000 
100 100-990 101-000 


For higher values of the line voltage, x, approaches the limit x) + 1 (see 
Figures 4 and 5), which means that at higher de voltages the amplitude of 


Figure 4 : Figure 5 
Maximum condenser voltage (*,) versus line Maximum condenser voltage (#,) versus line 
voltage (%)) for the exponential condenser voltage (%) forthe sinh-type condenser 


(for charging). (for charging). 
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oscillation is practically equal to the line voltage. At lower line voltages, 
however, the percentage increase of maximum voltage over the line voltage is 
greater. Comparing with the linear case in which x,/x) = 2 for all values of x9, 
we see that ¥,/x, is far less in the case of a non-linear condenser. This is a result 
of the strong nonlinearity. The capacity value rises so much with increasing 
voltage that the condenser cannot afford to have a large voltage across it, if the 
circuit equation [equation (2)] is to be satisfied. For exactly the same reason 
x,/X is greater in the sinh case than in the exponential case because sinh x/x 
is less than e* and hence the non-linearity of the sinh-type ee is less 
pronounced than that of the exponential one. 
From equation (10), differentiating the absolute value of the current 
through the exponential condenser, we get 


—— Slt ee AT ay eo al (x+1)=0 


atv = x, and 


a? |i | 
ae 


= ? [1 + e* (x % — x2 — 1+ x)]-*” €2@ (x 1)? (%,— 2)? 
+ 2 [+ 6 (x xy — x2 14 x)” 
x 6" ((% + 1) (%) — x) +(% —22%—1)], 


which is negative at x = Xp. Stnallarly, for the sinh-type condenser, we get from 
(15) 
d | } ‘ 
a ae [(cosh x — 1) + (% — x) sinh x]~+ (x, — x) cosh x = 0 


atx = x, and 


Ge | 
ali 


, 2 [(cosh x — 1) + (% — x) sinh x]~*? (x) — x)? cosh? x 


sh 2 [(cosh x — 1) + (x9 — x) sinh x]—1? 
x { (x 9 — *) sinh x — cosh x } 


which is also negative at x = x9. Hence for both the types of condenser | ¢ | 
has a maximum at «= x5, the line voltage. At x = x», therefore, 1 has a 
maximum if dx/dt is positive and has a minimum if dx/dt is negative. The 
situation is shown graphically in Figure 6 as well as in Figures 9 to 16. In 
Figure 6, both the points A and B correspond to % = x); however 7 has a 
maximum at A and has a minimum at B. 
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Figure 6 
Diagram showing typical charging waveforms. 
The values of the maximum current for different line voltages are tabulated 
in Table 2 and their relation is shown in Figure 7. 


Table 2 
Iman 
Xo 
Exponential sinh-Type 
1 1-000 0-737 
+ 12*837 5*+130 
10 445-20 104-94 
100 D188 x102% 3-64 x 1021 


Whereas in the linear circuit the maximum current is proportional to the 
line voltage, strong non-linearity is observed in the present case between the 


maximum current and the line voltage (Figure 7). 


charge:109 Imax 

— = nm 

uo OQ La) L—) 
discharge: 109 imax 


QoQ 


4 


Figure 7 
sinh; ——— exponential. 


Maximum current through the condenser (7,4) versus line voltage (x9) (for both types of condenser 
at charge and discharge). [Plotted on a log-log scale. ] 
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We have seen previously that from ¢ = 0, the voltage rises from zero and 
reaches a maximum at x = x,. Due to the convergence of the integral 


[ow ae, 


this maximum voltage is reached in a finite time, say (x,). After this time, 
the voltage must fall, since as previously proved, x has a maximum at %, 
and hence until the voltage falls to zero again, we have to use the negative 
value for dx/dt in equations (7b) and (13). In this range, the time equation 
becomes 


<M, ty * 
¢ \t/t oT ee oe ae 4 
= i @(x) dx a P(x) dx ef D(x) dx 


wy 


& (T(x) 
or 
a — I (x4) = I(x) — I(x). (26) 


Equation (26) is the equation of symmetry which expresses the fact that equal 
voltages are reached at equidistant times from the time of maximum voltage, 
i. e., that the waveform of x versus ¢ is symmetrical about the ordinate t = I(x). 
We can, therefore, plot x versus ¢ up to the maximum for x and then fold the 
curve about the ordinate ¢ = J(x,) to get the other symmetrical half. This 
procedure, when repeated at each zero of x, gives x as a function of time. 
Obviously x becomes periodic in time with period T = 2J(x,). This is to be 
expected since there is no dissipation of energy in the circuit. 

The values of the normalized time-period, found by numerical integration, 
are tabulated in Table 3 and plotted in Figure 8. 


Table 3 
£ 
%% 
Exponential sinh-Type 
il 18-286 Net) 
137-34 43-781 
10 4073-2 867-50 


| 100 4-096 x 1078 a 0285028 


Whereas in a linear circuit, the period is independent of the line voltage, 
in this case we observe that the former is strongly dependent on the latter as 
is shown in Figure 8. In Figure 8, for discharge, the values of the half-period 
are plotted instead of the full period, but in a logarithmic plot, this merely 
means a change of origin. 
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log: 


Figure 8 


sinh; ——-— exponential; a charge; b discharge. 


Period of oscillation versus line voltage (for both the types of condensers at charge and discharge). 
{Plotted on a log-log scale.] 


Discussion on the Discharging of the Condenser 


The convergence of the integrals (18) and (21) at x = x) and at x = —x, can 
be proved exactly similarly as we have done for the charging case. 

It has already been proved that x», the line voltage, is the maximum 
voltage across the condenser during discharge, i. e., that there is no ‘over- 
shooting’. In a similar fashion it can be shown that the minimum voltage 
across the condenser during discharge is —% 9. The result is quite analogous to 
the linear case. 

That the discharge is periodic can be proved exactly as in the charging 
case; but here the period is 4 J(x,), for to make one period complete the voltage 
has to fall from x, to zero and further to —%, and then again rise to %, in the 
reverse order. The periods of oscillation, T, at discharge have been calculated 
by numerical integration for some representative values of x) and have been 
tabulated in Table 4. A plot of this against x) has been shown in Figure 8. 


Table 4 
Hh 
X : ; wie 
Exponential sinh-Type 
1 14-000 9-517 
4 60-144 21-688 
10 1276-60 258-12 


bee 9-282 1022 | . 2-895x 1022 


It can be easily proved as in the case of charging that the voltage waveform 
at discharge is symmetrical about both the ordinates at t= 0 and at ¢= 21 (%9). 


_ ZAMP VIII/18 


274 ASOKENDU MOZUMDER ZAMP 


This waveform is further symmetrical about the axis of abscissae x = 0. This 
symmetry makes the average voltage in a full cycle zero which is essential in 
discharge. 


A study of equations (19) and (22) shows that the magnitude of the current 
is a maximum at x = 0 and has the value 


iy ev? Vx2 +1 — Xp 


for the exponential condenser and the value 


to Vx» sinh x) — cosh %» 


for the sinh-type condenser. The results are shown graphically in Figures 9 
to 16. 

In Figures 9, 10, 11 and 12 the waveforms of the voltage and the current 
for the exponential condenser for different values of the line voltage have been 
plotted by numerical evaluation of the integrals in equations (8), (10), (18) and 
(19). The waveform for a complete cycle is shown for charging and for a half- 
cycle for discharging. The close resemblance of the voltage waveforms to 
square waveforms is remarkable for higher line voltages. 

In Figures 13, 14, 15 and 16 the waveforms of the voltage and the current 
for the sinh-type condenser for different values of line voltage have been 
plotted by numerical evaluation of the integrals in equations (14), (15), (21) 
and (22). The waveform for a complete cycle is shown for charging and for a 
half-cycle for discharging. The close resemblance of the voltage waveforms to 
square waveforms is again remarkable for higher line voltages. 


Possible Applications 


Apart from their theoretical interest, the results obtained here may have 
some practical utility. For example it is seen from Figures 11, 12, 15 and 16 
that the voltage waveform approaches (both at charge and discharge) a square 
waveform for comparatively high line voltages, especially in the case of the 
exponential condenser. The effect is best observed at x) = 100, the highest line © 
voltage considered. This, therefore, indicates a simple way of generating square © 
waveforms of voltage by the use of non-linearity, provided that the voltage 
applied is not so high as to destroy the inherent non-linearity. 

Figures 10, 11 and 14 and to some extent Figures 9 and 15 show a current 
which approaches a triangular waveform. At x) = 100, however, in both cases, 
the departure from the triangular waveform is prominent. At this voltage the 
current waveform for charging is like a square wave. These may, therefore, be 
utilized to produce current waveforms of the types mentioned. 
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Figures 9 and 10 
a voltage at charge; b current at charge; c voltage at discharge; d current at discharge. 


Diagrams showing the voltage and current waveforms for the exponential condenser at charge 
and discharge for %) = 1 and 4 respectively. 
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Figures 11 and 12 
a voltage at charge; b current at charge; c voltage at discharge; d current at discharge. 


Diagrams showing the voltage and current waveforms for the exponential condenser at charge 
and discharge for %) = 10 and 100 respectively. 
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Figures 13 and 14 
a voltage at charge; 6 current at charge; c voltage at discharge; d current at discharge. 


Diagrams showing the voltage and current waveforms for the sinh-type condenser at charge and 
discharge for x) = 1 and 4 respectively. 
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Figures 15 and 16 
@ voltage at charge; 6 current at charge; ¢ voltage at discharge; d current at discharge. 


Diagrams showing the voltage and current waveforms for the sinh-type condenser at charge and 
discharge for %) = 10 and 100 respectively. 
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The dependence of the period on the line voltage can probably be used for 
producing frequency modulation, although this will be accompanied by ampli- 
tude modulation. The modulation will not be linear unless the applied voltage 
is small. 
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Zusammenfassung 


MacponaLp [1] hat gezeigt, dass die Kapazitaét eines Kondensators mit 
einem Halbleiter als Dielektrikum proportional zu 


sinha V 
aV 

variieren kann. V bedeutet die iiber dem Kondensator liegende Spannung, und 
« ist eine Konstante. ; 

Im folgenden wird eine Schaltung untersucht, bei der eine derartige Kapazitat 
in Serie zu einer linearen, verlustlosen Induktivitat liegt. Eine Methode wird an- 
gegeben, nach der Strom und Spannung als Funktion der Zeit fiir zwei Falle 
berechnet werden: 

1. Zur Zeit t = 0 wird die Kapazitat tiber die Induktivitat durch einen Span- 
nungss toss. aufgeladen. ; eas 
2. Die Kapazitat ist zur Zeit t= 0 aufgeladen und entladt sich liber die 
Induktivitat. | 

Die Schwingungsdauer, die Spannungsmaxima tiber dem Kondensator sowie 
die Maximalwerte des Stromes werden in Abhangigkeit von der angelegten 
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Gleichspannung untersucht. In diesem Zusammenhang wird auf die Moglichkeit 
hingewiesen, einfache Generatoren fiir rechteckformige und dreieckformige Span- 
nungen herzustellen. Die hier abgeleiteten Ergebnisse werden mit Rechnungen 
verglichen, die fiir einen Kondensator ausgefiihrt wurden, dessen Kapazitat wie 
exp(« V) von der Spannung abhangt. Derartige Kapazitaéten lassen sich mit 
einem Verstirker grosser Verstarkung, einer kapazitiven Riickkopplung und 
einem Netzwerk, dessen Ausgangsspannung exponentiell von der Eingangsspan- 
nung abhangt, realisieren [6]. 


(Received: November 26, 1956.) 


Uber naherungsweise Auflosung von Systemen 
homogener linearer Gleichungen 


Von ALEXANDER OSTROWSKI, Basel?) 


1. Die Lésung des linearen homogenen Systems 
mee Ai a (1) 


wo S eine ” x n-Matrix mit verschwindender Determinante ist und € als ein 
n-dimensionaler (Zeilen-) Vektor gesucht wird, kann nach den Vorschriften der 
Determinantentheorie ausgefiihrt werden, entweder durch die Aufstellung 
geeigneter Zeilen von Unterdeterminanten von S oder durch schrittweise Re- 
duktion von S auf eine kanonische Form, in der eine Reihe von Zeilen aus 
lauter Nullen besteht. 

Da aber in der numerischen Praxis sehr oft die Determinante von S nicht 
genau verschwindet, sondern nur «sehr klein» ist, ergeben sich dabei gewisse 
Schwierigkeiten und Unsicherheiten, die in zusammenfassenden Darstellungen 
iiber numerische Behandlung linearer Gleichungen eher nur fliichtig beriihrt 
werden. Es diirfte daher der im folgenden mitzuteilende Satz von Interesse 
sein, der eine allgemeine Methode zur Aufstellung solcher « angenaherter Lé- 
sungen » enthalt sowie eine Handhabe zur Beurteilung des gemachten Fehlers 
hefert. Wird namlich anstelle der singularen Matrix S die « gestérte» Matrix 
S+ A betrachtet, wo A = (a,,) eine m x n-Matrix mit gegen 0 strebenden 
Elementen 4,, ist, so wahle man einen festen n-dimensionalen Vektor 4 und 
betrachte die Gleichung 


(S +A) Bo aen (2) 


Verschwindet die Determinante von S + A nicht, so wird dann «im allge- 


1) Mathematische Anstalt der Universitat Basel. 
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meinen» | €| ins Unendliche gehen, zugleich wird der normierte Vektor &/ (Ee 
sich von einem geeigneten Lésungsvektor von (1) um 


O (Vd | a, v | ‘ 
My? 


unterscheiden. 


2. Genauer lautet unser Ergebnis wie folgt: 


Satz 1. Es set S eine n x n-Matrix vom Rang n—k,1<k<n—1. Es 
set A = (a,,,) eine n x n-Matrix mit variablen Elementen a,,,, die gegen O streben, 
derart, dass dabet die Determinante von S +- A von 0 verschieden ist. Es sei n + 0 
ein konstanter V ektor, fiir den die Gleichung 


So = 9! (3) 


keine Losung hat. Dann gibt es zu einem Liésungsvektor E von (2) einen (variablen) 
Lésungsvektor x von (1), derart, dass 


Spas Jy lA" ¥ 


My? 


Orn! Srp 


gilt2), 


3. Beweis. Bekanntlich lassen sich zwei nicht singulaére konstante n x n- 
Matrizen B,C derart finden, dass 
Gg. O 
BSC =S,= (5) 
O En-x 


gilt, wo O, eine aus lauter Nullen bestehende k x k-Matrix und E,,_, die Ein- 
heitsmatrix der Ordnung » — k ist, wahrend die beiden Matrizen O aus lauter 
Nullen bestehen. Setzen wir dann 


BAC = Ag = (a) » (6) 
COPE Cyt Me Sey pe Bay = Nig : (7) 


so gehen die Gleichungen (1) und (2) tiber in 
So ah oat (1) 


(Sy +A) Bam. (2*) 


2) Die Bedingung, die im obigen Satz dem Vektor 7 auferlegt wurde, lasst sich bekanntlich 
auch so formulieren, dass der Vektor 7 nicht auf der gesamten Schar der linksseitigen Nullvektoren 
yon S senkrecht stehen darf. In der numerischen Praxis lasst sich dieses Kriterium allerdings nur 
sehr selten verwenden, so dass man eventuell fiir 7 nacheinander m linear unabhangige Vektoren 
einsetzen wird, zum Beispiel die » Koordinaten-Einheitsvektoren. 
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Setzen wir dann 


No = (Vis s>d Yn) so (Pay reke oe daalt (8) 
so gilt nach Voraussetzung des Satzes fiir ein geeignetes x zwischen 1 und k: 
yO, ol SW Sek (9) 

4. Wir setzen ferner 
(Kise exe Oy «22 Ose, (Ope On ere eee (10) 


so dass =o + 9 ist. Offenbar ist Syo’ = 0, Sy &, = 0. Um Ay &, geeignet dar- 
zustellen, fithren wir allgemein den der v-ten Kolonne von A, entsprechenden 
Vektor C, ein, so dass 


Gi ce (gO) ie, FON s lack ea teeeages 


7 NY 


ist. Dann kénnen wir (2*) in der Form 
(So + Ap) £9 = D/ %, C, + 0” = No (11) 
y=] 


schreiben. Vergleichen wir rechts und links in (11) die x-te Komponente, so 
folgt 


n 
= (0) 
LS o> any Xy > 
(Dil, 


da x <k ist, und daher 


lel Sol Vee. (12) 


5. Setzen wir daher fiir die Norm | Ay |, von Ay: 


a = |Ayla= |/2y bane? (13) 


so folgt wegen (9) 


1 a 
ey eae |&| > 00. (14) 
Dividieren wir nun die rechtsseitige Relation in (11) durch | & | , so ergibt sich 
n C, erg est : 
Tet ot Tet = Tar? - 


und daraus, da or |C,|* =a? ist, nach der Cauchy-Schwarzschen Ungleichung 
v=1 
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wegen (14) Ht he 
1¢l. io <a ‘No | 
a ot = Teh te Se ler 1). (16) 
Damit gilt nunmehr 
f9 =o + @—o+ O(\Eo| a) . (17) 
6. Aus (17) folgt weiter wegen (7) &'=Co'+Co', und daher, wenn 
Co’ =| &| x’ gesetzt wird — unter x eine geeignete Lésung von (1) verstanden -, 
c= | ela + C oo tex’ + 0 ((E5| a) . (18) 


Aus (7) folgt aber 
[Fol = O(/ |) , 


so dass aus (18) durch Division mit | €| endlich 


" 


Rs zw + O(a) (19) 


[Ale= 7D lapel? 
MY 
und bilden in analoger Weise | B |, und | C |,, so folgt bekanntlich®) aus (6) 
a = |Aole =|Bl. [Ale (Cle 


wre 


folgt. Setzen wir nun 


und daher fiir ein konstantes c = | B |, | C |o: 
aes clAl,. 


Tragen wir dies in (19) ein, so ergibt sich die Behauptung (4), und Satz 1 ist 
bewiesen. 


7. Der Vektor z in (4) kann « wesentlich variabel» sein, das heisst braucht 
auch nach geeigneter Normierung nicht gegen einen konstanten Lésungsvektor 
von (1) zu konvergieren. Man betrachte zum Beispiel fir u > 0, v > 0 die 
Matrizen 


0 0 0 u 0 0 
Sieh <0 0 hee (Sas Ave} pod ~ 0 
O01 Oa) vi 


Fiir 7 = (1, 1, 1) ergibt sich € = (1/u, 1/v, 1), 


i 1 
z SS 
u? v il 
Jrtwe Jitee e i+ dete 
3) Vergleiche zum Beispiel A. OSTROWSKI, Uber Normen von Matrizen, Math. Z. 63, pp. 3 und 
13 (1955). 
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Streben nun u und v derart nach 0, dass fiir geeignete Teilfolgen sowohl u/v > 0 
als auch v/a > 0 gilt, so ergeben sich in der Grenze als zugehérige Haufungs- 
vektoren von (20) sowohl (1, 0, 0) als auch (0, 1, 0). 


8. In dem obigen Beispiel hangt A analytisch von zwei gegen 0 strebenden 
Parametern ab. Handelt es sich dagegen nur um einen analytisch eingehenden 
Parameter, so konvergiert a nach geeigneter Normierung in der Tat gegen 
einen konstanten Lisungsvektor von (1) von der Lange 1. Genauer: 


Satz 2. Unter den V oraussetzungen von Satz 1 seien die a,, analytische Funk- 
tionen eines Parameters t, reguldr in t = 0. Setzt man dann 


es (21) 


so gilt fiir den Vektor x in (4) mit einem geeigneten natiirlichen 4 
fat=t+O0(t), (22) 
wot ein konstanter Losungsvektor von (1) von der Lange 1 ist. 


9. Beweis. Wird (2) aufgelést, so erhalt man fiir den Vektor €, nach wach- 
senden Potenzen von ¢ entwickelt, 


&=14(U+0()), (23) 


wo U ein von ¢ unabhangiger konstanter Vektor ist. Der Exponent 7 ist hier 
positiv, da | € | ins Unendliche strebt. Daraus folgt weiter 


= «i (r+ 0(i)), (24) 


wo U/|U|=t gesetzt ist. Andererseits folgt aus unseren Annahmen, dass 
ayy = Olt), | A |, = O(é) ist; und (22) folgt nunmehr aus (24) und (4). 


10. In der Praxis ist wohl am wichtigsten der Fall, wo A eine skalare 
Matrix ¢ E,, und ¢ der « Fehler» des in Betracht kommenden Eigenwerts ist. 


Dann ist | A |, = tV/n, so dass der Fehler bei der Bestimmung des zugehérigen 
Eigenvektors von der gleichen Gréssenordnung ist wie derjenige des Eigenwerts. 
Andererseits ist in diesem Falle die natiirliche Zahl 7 der Formel (22) gleich 1, 
sofern S nur lineare Elementarteiler hat, also insbesondere reell-symmetrisch 
oder allgemeiner hermitisch ist. Es folgt dies daraus, dass, wenn der charakte- 
ristischen Wurzel 0 von S nur lineare Elementarteiler entsprechen, die Deter- 
minante von S+#E fiir t= 0 eine Nullstelle genau der k-ten Ordnung hat, 
wahrend ihre Unterdeterminanten der ersten Ordnung durch ?¢*-1, nicht aber 
durch ¢” teilbar sind. 
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Ist insbesondere S reell-symmetrisch und ¢ reell, so folgt 


= tT sgni + Off). (25) 


re) Ore 


| 
Dies liefert in der Praxis die Méglichkeit, die Konvergenz gegen den Vektor t 
zu erhalten, indem man durch Vorzeichenanderung die erste nicht merklich 
gegen 0 konvergierende Komponente von &/|&| durch eventuelle Multiplikation 
mit —1 positiv macht. 
Summary 


Ifa singular n xX m-matrix S is approximated by S + A, an approximation 
toa non-trivial solution of S¢’ = 0 is obtained from a solution £ of (S + A) & = 7’ 
in forming &/| €|. The vector 7 is here a conveniently chosen constant vector, 


and the error is, if d = (a 
i 
0 (VX 4 HY 


BY 


lt P ~ 


(Eingegangen: 29. Dezember 1956.) 


Spannungsfeld der auf den Rand einer elastisch anisotropen 
Halbebene wirkenden Tangentialkraft 


Von Hans H. STADELMAIER, Raleigh, North Carolina, USA14) 


Einleitung 


Das Spannungsfeld einer Normalkraft, die auf eine elastisch anisotrope 
Halbebene einwirkt, wurde in einer kiirzlich erschienenen Arbeit?) beschrieben. 
Die darin gewonnene Lésung gilt fiir beliebige kristallographische Systeme, 
mit der (aus Griinden der mathematischen Vereinfachung eingefiihrten) Ein- 
schrankung, dass die x- oder y-Achse eine geradzahlige Symmetrieachse sein 
muss. Die Ubersichtlichkeit des dort gewonnenen Ergebnisses regte zu der 
nachfolgenden kleinen Untersuchung an, die den Einfluss einer parallel zum 
Rande verlaufenden Einzellast beschreibt und schliesslich mit dem friiheren 
Ergebnis zu einer Randlast beliebiger Richtung kombiniert werden kann. 


Problem und Loésung 


Wie man aus dem Lehrbuch von VorerT?) entnehmen kann, lautet bei gerad- 
zahliger Symmetrie der x- oder y-Achse die Differentialgleichung der Span- 


1) Department of Engineering Research, North Carolina State College. 
2) H.H. Stapevmarer, ZAMP 7, 393 (1956). 
3) W. Vorct, Lehrbuch der Kristallphysik (Leipzig und Berlin 1910). 
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nungsfunktion at at of 
Pp YP ee 
pp a ea eye Ore (1) 
wobei 
a= Sia + S60/2 und b ae Ste 


$14 Si4 
und die s, Elastizitatskoeffizienten in der Voigtschen Bezeichnungsweise sind. 
Der Strich deutet an, dass die Elastizitatskoeffizienten von der Orientierung 
zwischen (x, y,z) und den kristallographischen Achsen abhangen. Die aniso- 
trope elastische Platte erfiille die Halbebene y 2 0. Die Randbedingungen 
lauten jetzt 


Try (x 0, y= 0) =0,. Tangentialkraft, P in«(0,0).% 0, @,0) =Oeaat 


(1) und (2) werden erfillt von der nachfolgenden Spannungsfunktion, die in 
zwei Abschnitte zerfallt : 


Falls a* > 0, 
P 
e ap—a) al? y arctg ?» igs + In (x? + p? y2) + xInp 
(3a) 
= qy arctg-4* t+ 5 ae n(x? + g? y?) — x Ing] 
wobel 
p>=a—Var—b, ge=atVar—bd 
bzw. bei a2 < 5b, 
rae le vY ry 
i ates {(#—s 9) arctg "9 spe eae ale (3b) 
= ee 2 2 
(x + sy) arctg—— 3 Infr? 9? + (x +s »)*I} 


mit 
Tris S (a+Vo), 3? = . (—a + Vo) ? 


Fir f= q=1 bzw.7 = 1, s = 0 geht (3 a, b) tiber in den bekannten Ausdruck 
fiir den isotropen Fall 


P 
y= — yaretg. (3c) 


Aus der Spannungsfunktion » folgen durch Differentiation die Spannungen 


heey Op ay? 

1 oat oe PTD Gee pe yy GEE GTO’ 
nie, Ope ar evi 

ee Oy2 og (b + 9) (v2 + p? y2) (42 + g? yp)’ (4a) 
ee ees wy 

a dx dy ~ me +9 Carpe y?) (#2 + gq? y?) 
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bzw. 

2.P ry? 

rd ot aaa ae ae ey Mn a 
a yt te — sy) Pye + (et sy) h]’ 

ee 2P , x 

Ta” PE NE et OT Be 

2P p2 ay 

Tay =e aie tia. 


mE ey Hat NT 
Fir = q= 1 bzw. r= 1, s = 0 erhalt man wieder den bekannten isotropen 
Fall 


eer ee tT oe Si ee ewer ; 
i RL (x + ape) 2 p & It (a2 a5 apehie / vy y?) ge ( c) 


Transformation auf Polarkoordinaten 


Die Beziehungen y= Rsin#, x = Rcos# fiihren unter Beriicksichtigung 
der Transformationseigenschaften der Spannungen zu den Ausdriicken 


- BG Ja) So's ) cos 1 
a q sint # (ctg?& + p2) (ctg2o+ gq?) R’ Sa 
Og =0, Tra = 0 
bzw. 
ye ae : cosé Ne oat 
ia a sin! # [v2 + (s — ctg#)?] [72 4+ (s+ ctg#)2] R’ (5b) 


Og= 90, Trep=O0- 


Dies fiihrt fir p = g = 1 bzw. 7 = 1, s = 0 zum isotropen Fall 


PY Cos: 


I ie (5 o 


Op = 
Wie im Falle der Normalkraft hat man eine reine Radialspannung, welche fur 
die beiden friiher?) besprochenen Falle in Figur 1 dargestellt ist. Dazu muss 
man die s/, fiir die interessierenden Orientierungen in «-Eisen berechnen und 
daraus die ~, g bzw. 7, s gewinnen. Die Orientierungen sind (unter Angabe der 
Millerschen Indizes derjenigen Kristallachsen, die mit x, y und z zusammen- 
fallen) wie folgt: 


PalGea) x= <0), y—dlo), 2= 007), 
Fall. (5b) ~—100)., .y =—(10), 2= 110). 


Dafiir wird im ersten Fall # = 0,47, g = 2,12 und im zweiten Fall 7 = 0,47, 
s = 0,54. Mit diesen Werten sind die og berechnet. 
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i 60° 


Figur 1 


Polardiagramm der Spannungen nach Gleichungen (5) fiir @ < 90°. Kurven (a), (b) und (c) ent- 
sprechen Fallen a, b und c des Textes. 


Randlast in Schrégrichtung 


Fir eine in beliebiger Richtung unter dem Winkel « zur x-Achse (Figur 2) 
wirkende Randlast K gewinnt man das Spannungsfeld durch Superposition 


K Ksine 
Figur 2 
Lage der Winkel # und a. 


des Feldes der Tangentialkraft K cos« auf das der Normalkraft K sing. Fiir 
eine Normallast P hatten wir frither?) gefunden, dass 


il 1 


12 
Opie a p q (p cc q) sin? & (ctg?d ag p?) (ctg? 8 at q?) : R? (6a) 
bzw | : 
Pye) 
Og ae Troe St) : = (6b) 


sin? } [v2 + (s — ctg 8)?] [v2 + (s+ ctg#)7] R- 
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Damit ergibt sich fiir die Randlast K 


15 Pq sind sing + cos # cos « 1 
oO = —_ —— —- - — — ~~ ¢—— 
R BL (+4) sint } (ctg?d + p?) (ctg?9 +942) R’ (7a) 
bzw. 
ak y2 <2 j j . o Vets 
grenk ‘ kaa SSA B08 feos ¢ COSa tal te (7b) 
I sint # [y? + (s — ctg @)?] [v2 + (s+ ctgd)?] R 


Fir Isotropie (6 = g = 1 bzw. r= 1, s = 0) erhalt man wieder den bekannten 
Fall 


DK. ob 

Ce eee Aa (70) 
Betrachtet man den gesamten Nenner in (7a) bzw. (7b) als verzerrten Radius, 
so bemerkt man, dass die Superposition der Spannungen zu Ausdriicken fiihrt, 
die im Gegensatz zum isotropen Fall nicht mehr einfach aus der Differenz der 
Winkel « und # hervorgehen. Dies ist eine Folge der unterschiedlichen Koeffi- 
zienten in Gleichungen (5) und (6). (Die Koeffizienten miissen die angegebene 
Form haben, damit die Randbedingungen erfiillt werden.) 


Diskussion 


Wenn man bedenkt,.dass fiir Fall a der Figur 1 die Richtungen maximalen 
Elastizitatsmoduls parallel x und y verlaufen, fiir Fall 6 dagegen bei # = 55° 
und 125° legen, so findet man wieder das Prinzip bestatigt, dass die Spannung 
in den elastisch harten Richtungen verstarkt wird. Insbesondere findet man 
fiir die schragliegende harte Richtung (Fall b) die maximale Spannung wie- 
derum zwischen der hartesten Richtung und der Richtung der Randkraft. 

Im Problem der Normalkraft wurden fiir Fall b die Héchstwerte der Span- 
nungen bei #= arcctg(+s) gefunden, also in der Nahe der Minima des 
«reduzierten Radius»*) 


n= wey sin4d [72 + (s — ctg)?] [v2 + (s+ ctg#)?], 
wahrend man mit der Tangentialkraft, trotz gleichen Wertes des reduzierten 
Radius, die Maximalspannungen starker zu der x-Achse geneigt findet. Dies 
erklart sich daraus, dass die aus 0op/00=0 und ctg?—-+s abgeleitete 
Bedingun 
mer = (377+ 8s?— 4) =0, 


deren linke Seite ein Mass des Anstiegs an den Stellen # = arcctg(+s) ist, 
von den Zahlenwerten fiir y und s mit guter Naherung erfiillt wird. Dagegen 


4) Nur ausnahmsweise wird in einem dem Fall b entsprechenden Beispiel das Minimum des 
reduzierten Radius genau in der durch die Kristallsymmetrie ausgezeichneten <111>-Richtung 
 liegen, da neben der Symmetrie auch die Zahlenwerte der elastischen Koeffizienten hierfiir mitbe- 
stimmend sind. 


ZAMP VIII/19 
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findet man bei der Randbedingung der Tangentialkraft einen entsprechenden 
Ausdruck, der fiir arcctg(+ s), das heisst # < 90°, negativ und fir arcctg(—s), 
das heisst > 90°, positiv ist. Daraus folgt, dass bei der Randbedingung der 
Tangentialkraft die Richtung maximaler Spannung grundsatzlich zwischen 
der Richtung der Randlast und der Richtung des kleinsten reduzierten Radius 
liegen muss. 

Diese Arbeit wurde geférdert vom Office of Ordnance Research, US Army, 
wofiir an dieser Stelle gedankt sei. 


Summary 


The stress distribution obtained by solving the two-dimensional problem in an 
anisotropic medium, with boundary conditions of a concentrated tangential load 
acting on the boundary of a semi-infinite plate, is purely radial. The solution is 
given in closed form and is combined with the solution for a concentrated normal 
load to solve the problem of an inclined force acting on the boundary. 


(Eingegangen: 2. April 1957.) 


Toroidal Oscillations of a Spherical Mass of Viscous 
Conducting Fluid in a Uniform Magnetic Field 


By KeItH STEWARTSON, Bristol, England?) 


1. Introduction 


The torsional oscillations of a highly conducting sphere in a uniform 
magnetic field have been considered by PLUMPTON and FERRARO [4]?). They 
showed that if the conductivity o were assumed to be infinite in the governing 
differential equation but the boundary condition on the surface of the sphere 
appropriate to a finite conductivity retained, then a continuous spectrum of 
eigenvalues could be obtained. Corresponding to each eigenvalue there was an 
oscillation of one shell ot the sphere, formed by rotating a certain line of force 
about the axis of symmetry. The continuous spectrum is a novel feature in such 
problems and it was suggested by CowLtnc [1] that this degeneracy may be ~ 
removed if a rotating sphere is considered. 

In this paper we show that if the fluid is not quite a perfect conductor 
then there is only a doubly infinite discrete set of eigenvalues and the cor- 
responding oscillations are confined to the immediate vicinity of the axis of 
symmetry. If now o > oo the set remains discrete but becomes singly infinite 


1) University of Bristol, Department of Mathematics. 
2) Numbers in brackets refer to References, page 297. 
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and the corresponding oscillations occur only on the axis of symmetry. Any- 
where else in the sphere the fluid is at rest. The effect of a small viscosity is 
also considered and found to be negligible. 


2. General Considerations 


It is supposed that the sphere has centre Q, is of radius a and the magnetic 
field imposed is uniform, of strength H, and acts in the direction O z. If P is 
any point let R denote the distance of P from O, ry its distance from O z and 0 
the angle between O z and O P. Further if v,h are the velocity and magnetic 
vectors at P in the disturbed state MAXWELL’s and EULER’s equations reduce to 


oh 1 
oe ate curlcurlh = curl(u , A) , (2.1) 


, er 


=, + (v-grad)v = —— gradp + soe curllh j»h—vcurleurlv, (2.2) 


ole 


divu = divh=0 (2.3) 


where o is the conductivity in electro-magnetic units, @ the (constant) density, 
p the pressure, v the kinematic viscosity and where all body forces save magneto- 
hydrodynamic ones have been neglected. 

We assume that in the disturbed state all perturbed components oscillate 
with the same period 27/« and are small. For toroidal oscillations therefore, 
symmetrical about O z, we may write, in cylindrical polar coordinates, 


v= (0,rQé,0), h=(0,rhei*, H) (2.4) 


where Q, h are small functions of 7 and z only. The governing equations then 
reduce to ~p = const, 


: 02 it 07h st ey? 07h 
tah—= H + aaa (Gp yor st (2.5) 
and e ie 
H Oh 022 3 0) O 
; : : Bie tip Ve 6 
sions eae? ee. ny a) es) 


There are two boundary conditions at the surface of the sphere R =a. 
First h must be equal to the external irrotational magnetic field induced by 
the motion, and since the potential of this field must be a function of 7 and z 


nl ; 
at h=0 when R=a, (2.7) 


Second the stress across a tangent plane to the surface of the sphere at the 
sphere must be continuous and since it is constant and normal to the plane 
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ifR>a 
oN 2 Dagon alee (2.8) 


OR 
The problem is to determine which values of « lead to nonzero values of Q 
and h satisfying (2. 5)—(2.8). This problem is not easy in general, but greatest 
interest centres on the special case when a is large and y is small. However, if 
we formally let o > oo and y + 0 in (2.5) and (2.6) contradictions are obtained. 
For then the equations reduce to 


hehe wendune Ge ae 


Oz 4n0 Oz (29) 


with solution 


(4 0) Q=71 A(r)sin(Bz+e), h=A(r) cos(Bz + @) (2.10) 


where A and « are functions of 7 only, and 47 9 a? = H? f?. 

First from (2.5) and (2.6) we see that 02/07 and Oh/dr must vanish at 
vy = 0 for allo, v, and z. Hence if 2 and / are not identically equal to zero, they 
cannot both vanish everywhere on the axis of symmetry. In the limit ¢ > oo, 
vy > 0 therefore this must also be true implying that 


Al} 0" (2.11) 


Second, since A = 0 when R= a 
AQ) = OMT 7, dO oF + , and cos[f (a2 —72)\"+<6/=0 (212) 


where 7) < a is an arbitrary constant. 
This is PLuMpTON and FERRARO’S [4] result that the shells y = const. oscil- 
late independently and we note that it contradicts (2.11) unless 7, = 0. 
Third, since 0Q2/0R = 0 when R=a 
e=0 os, cos(B a+ «) =0 
and (2713) 
vy A'(r) cos(B z+ €) = B A(r) zsin(B z+ &) 


where z = (a? — 7?)"?, which contradicts (2.12). 

One reason for these two contradictions is that proceeding to the limit 
of zero viscosity or infinite conductivity or both, the order of the governing 
equation has been depressed by two so that two independent solutions have 
been lost. Although in general they will be negligible when either of » and o} 
are small, they will become important near the surface of the sphere since the 
other two solutions cannot satisfy both boundary conditions. A boundary layer 
will therefore be formed near R = a in which the operator 0/0R will be large to 
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counteract the numerical factors » or o~! in (2.5) and (2.6). The operator 0/00 
on the other hand can only be large in a small region of the surface for other- 
wise the oscillation would be too complicated to be of interest. 

Since if 0/OR is large so is 0/0z, both i ah and i «Q may be neglected in 
(2.5) and (2.6) which then reduce to 


el, 


eee OE: 1 a2h dh 
H cos6 spt = Hcos0 5 +42 0% >7: 


bpee SE . = 0. » (2.14) 
The solution of these equations which is rapidly varying when a — R is small 
positive is 


h = —4.2 (00 )"2.Q = C(6) exp \H (R — a) () cos 0} =h, (2.15) 
where C(@) is an arbitrary function of 6. The other rapidly varying solution 
corresponds to a boundary layer just outside the sphere. 

This argument applies equally well whether y or o~! or both are supposed 
small, provided only that neither is large. We shall now show that if v is small 
the boundary condition 022/0R = 0 is accounted for by a solution of the form 
of (2.15) so that in the limiting case » > 0 it may be neglected. For if y is'small 
the slowly varying solutions are found by expanding in powers of v. Let us 
call the leading terms, which satisfy (3.1) below, 4, and Q,. Then the complete 
solution is ' 


h=ht+h, and 4n72=422,—h (9 0%)7"” (2.16) 
and from the boundary conditions 


bits re cont C(O) = 0 (2cLn 


hy + C(8) = 9, ‘OR 4270 
when R = a. Hence 
H cos ln +4000 a0 == () (2.18) 


when R = a, and on proceeding to the limit vy > 0 it follows that hy = 0. 

According to this argument therefore the boundary condition on / must be 
satisfied by those solutions of (2.5) and (2.6) which neglect » as a first approxi- 
mation while the condition on 2 introduces a boundary layer in which 0Q/0R 
falls rapidly to zero on the surface but in which the changes in Q and h are 
only O(v g/o)"”. A rapid change in h in a boundary layer near the surface 
would of necessity result in a much more violent change in 02/0R so that if 
a boundary condition is to be lost as y > 0, it must be that on 0Q/OR. 

The solution in (2.13) is therefore incorrect, thus resolving the contra- 
diction between it and (2.12). There remains the contradiction between (2.11) 
and (2.12). This may be resolved if either there is a boundary layer in the 
neighbourhood of r = 0 so that the condition A(0) + 0 may be neglected and 
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which leads immediately to PLuMpTON and FERRARO’S result [4], or by requiring 
that 7, = 0. We shall show in the next section that the second alternative is 


correct. 
3. Eigenvalues for an Inviscid Conducting Sphere 


In the previous section we showed that if » is small the fluid could be 
treated as inviscid but having a finite conductivity in determining the eigen- 
values. Then h/ satisfies 

2 ) ) 2 \ 
P a He eS bee (2 ha bag ee ee 3.1) 


oz" \4n0 a ae0: 


and 4 = 0 when R= a. 
Introduce quasi-elliptic co-ordinates §, 7 where 


Y= (—)" cosh & cosy , z=b (; cit + a)" sinh & siny (3.2) 


bis aconstant and z = e*’”, Then the sphere R = a becomes & = & where 


: = tH@ i aH 
A sinh €, = ( A Iv ", provided that b ( . ay [ ea and 
o | 420 4nof 
07h < fi Oh; 2 2 
Fy y+ 3 tanhé SE — 3 tan Vise eg (cosh 2 € — cos2 n) me 


Since h=0 when &= &, the eigenvalues may be found by separating the 
variables. Write . 


h= X(n) Y(é) (3.5) 
when 
d2X cas a2 b? ; 
dat ~ 3tann G+ (c— => cos2q) X= 0 (3.6) 
and 
azV 
Gandia nhé _ ~+(%5 ~ cosh2 & — c) Y a0), (3.7) 


c being at present an arbitrary constant. 
The boundary condition now becomes 


Y=0 ‘when €= 4). 


Further if Y(0) + 0 then X must be an even function of 7 and Y’(0) = 0 so 
that h and 0h/0z may be continuous at z = 0. If Y(0) = 0 then X must be an 
odd function of 4 so that 0h/0z may be continuous at z = 0. Finally X must be 
bounded when 7 = +2/2. 

First let us consider the equation for X assuming to begin with that « b and 
c are real. If either X(0) or X’(0) > 0, X must increase steadily in 0 <4 < nq 
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where «? 6? cos 2 9 = ¢, oscillating thereafter but remaining finite as 7 > 2/2. 
In particular if X only oscillates a finite number of times then as « b +0, 
No > 2/2. In that case the conditions imposed on X at 7 = 0 imply that, if 
X(z/2) = 1, whichever of X(0) and X’(0) does not vanish must be exponentially 
small. An approximation to the acceptable values of c and the corresponding 


forms of X as « b oo will now be found. First of all consider values of y near 
m2: Let 


"3-5 0) 
when ‘ } . ‘ 
Gh saths(ax-d4) 09 


for finite x where 8 a b B = 2c + a? b?. 
The boundary conditions are that X = 1 when x = 0 and X is exponentially 
small when 7 is small and therefore when x is large. It is now easy to verify 
that this is only possible if 
pan+1+0(-,) (3.10) 
where ” is an integer. 
The solution is then 
. (—)? n! x* 
(n—r)!r! (r+ 1)! 


See a 
r=0 


1 
+ a ) (3.11) 
and is valid so long as « = o0(« 0). 
On the other hand, so long as 
<=> @b)-¥ 
is a solution of (3.6) in descending powers of « b may be obtained by JEFFREY’S 
method [2]. Write 


X = F(n) expa b (siny — 1) 
and substitute into (3.6) assuming that F is a slowly varying function of 7. 
The first approximation gives 


cos?” 1 
_ These two solutions (3.11) and (3.12) overlap and may be combined into one on 
writing 
x = 2ab(1—siny) 
- instead of (3.8), viz. ; 


xX —— 


gan+2 exp b (sin 4 = 10) os (—)" n! (2 a b (1 — sin 7)I” . (G:15) 


(1 + sin )2”+2 = (n — r)!r! (x + 1)! 
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Since this solution is exponentially small when 7 = 0 it follows that only 
very small multiples of that solution of (3.6) which is exponentially large need 
be added to it to make either X(0) or X‘(0) vanish. Since an addition of that 
kind is made by altering the value of c we conclude that to the order of the 
approximation there are both odd and even solutions of (3.6) bounded at 
n =7/2 when 

canto 
Ce meee 4(n+1)a+O(1). (3.14) 


A similar result is known for Mathieu functions, and the reader is referred to 
(3] for details. 

Second let us consider the equation for Y where ¢ is given by (3.14) and 
either Y(0) or Y’(0) = 0. Since only o determines the order of magnitudes of 
the parameters £, = O(b-!) we may expand (3.7) in powers of &, retaining only 
the leading terms, and obtaining 


OY + [a2 B® — 4 ab (m+ 1] Y = O(Y a2 62 E4), (3.15) 
whence 
Y= Asin(§[a?b?—4ab (m+ 1) +e), (3.16) 


where A and « are constants. 
Hence from the boundary conditions 


Ey [a2 b2 —4ab (n + 1)]}'? = ara + 0(é) (3.17) 


where m is an integer. 

Although (3.17) has been obtained on the assumption that « b is real we 
may now appeal to the theory of analytic continuation to extend it to the 
values given in (3.3) provided only that we avoid « = 0. It then follows that 
when o is large 


mH xa (ee ae | m2 H2 
2a 2a \|42ea%c? 


i. (eae o(-) (3.18) 


where m and are integers. 

For finite values of m and n the oscillation is confined to the immediate _ 
neighbourhood of the axis of symmetry being negligible at distances O(o—”) . 
from it. Hence in the limit ¢ > oo the disturbance is confined to the axis of 
symmetry so that the criterion (2.11) is of equal importance to the condition 
h=0 when R= a. 

When o is large and x is finite (3.18) shows that the oscillations have a 
damping factor O(o—). It is possible to obtain solutions in which the oscillation 
is not confined to the neighbourhood of the axes but only if 2 = O(c") from 
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(3.14) and (3.3). In these cases however the time of decay is O(a/m H) which 
is of the same order as the period of oscillation. 
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Zusammenfassung 


Es werden die toroidalen Schwingungen einer fliissigen Kugel in einem kon- 
stanten Magnetfeld betrachtet. Nach [2] ist bekannt, dass das Spektrum der 
Eigenwerte fiir vollkommen leitende Fliissigkeit kontinuierlich ist. Hier wird 
aufgezeigt, dass fiir endliche Leitfahigkeit ein diskretes Spektrum entsteht und 
_dass der diskrete Charakter erhalten bleibt, wenn die Leitfahigkeit nach Unend- 
lich strebt. In diesem Grenzfall konzentriert sich die Schwingung auf die Sym- 
metrieachse. 


(Received: December 22, 1956.) 
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Résolution systématique de problémes aux limites linéaires 


Par Paut Bureat, Neuchatel 


1. Introduction 


Dans un précédent travail [1}+), j’ai indiqué la transformation fonctionnelle 
appropriée a la résolution, dans un intervalle fini, de problémes constitués par 
une équation différentielle du m-iéme ordre et m conditions aux limites générales, 
linéairement indépendantes, et donné la solution sous la forme d’une série de 
fonctions propres. 

Rappelons quelques résultats relatifs au 2° ordre. 

La méthode appliquée au probléme 


y" +r y= f(x); oy (0) + 4 90) = 0, By (2) + By y'(!) = B 
({(*), fonction généralement continue; y, les « et les 8 sont des constantes réelles ; 
(2 + o3) (B23 + Bi) #0; a, %, 50; Bs By = 0} conduita 


= a OK a 
= —— | - — A. SIN @,, ¥ — Oo Wz COS WM, *% 1 
y Den mer os (4 — 0 B) — Ag] (1 Sin oy ¥ — ay C08@% 2) (1) 


1) Les chiffres entre crochets renvoient a la Bibliographie, page 302. 
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1 2 wr— at . ; 
on Ny = (ae + a) wh) SF Ma Ok Pein? My 1 — a & Sin? w, 1 
ed, 4 Ww, 
= (a2 + a8 w2) l bil (a2 By — %1 Bs) (%. By ver at Bs) 
= 2 GGA Nae Same eee eAn bs 
te ot me 2 (BR + Bi oF) 
a, SIN Wp — Hy Wy COSM,1 _ %, COS oy, | + % Oy, SiN Or L 
a Ba Ox Bs 
Oe 
1 
1 : 
A= N, ‘ T(&) (a, sin wy, & — %, Wy COS Wx &) dé; 
IC 


0 
les w, sont les racines positives de l’équation 
k Pp 


oe Ba ois papal sc ate eae 7 (2) 
oo Bs — a By a Bs — a By wo 


Quand a, = Bs = 0, a = B, = 1, la série écrite plus haut doit étre précédée 
du terme , 
1 }a—f68 it ¥ 
po + fue al. 
0 


Cas particuliers 


cotg w 1 


AO) = 10.0). y= )/ - ea ( 1) 6] - 4, | sink 7 95, 
bien) (3) 
I 
4, = [ a) sink T Eas 
0 


Xx + ot a IE oe A, cos he x: 
f 1 
Ay = SG GE; AG == [He) cosk7-§d5 (k + 0) 
: 0 
y'(0) = 0, ¥(1) = B: y= — Le 
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[ear 


4 |<; ((24 4) a — (-1)4 2 f) —A,| sin(# — 7) + #3 


a! 


Remarque. Les séries ont pour somme y en tout point de l’intervalle ouvert 
(0, 7). Si a By + 0, elles représentent la solution sur (0, /). 


2. Equations aux dérivées partielles 


Soit a résoudre un probléme aux limites 4 deux variables indépendantes, 
Vune d’elles variant sur un segment, l’autre étant le temps ou, plus généralement, 
une variable parcourant l’intervalle (0, + oo)?). 

Habituellement, on commence par appliquer la transformation de LapLaceE; 
on obtient un probléme aux limites a une seule variable indépendante que l’on 
résout. La solution constitue la transformée de LapLtace du probléme proposé; 
on termine par une inversion. 

Si ce procédé ne convient pas, soit que les calculs conduisent, lors de l’inver- 
sion, a des contours d’intégration compliqués [2], soit que les conditions initiales 
requises pour l’emploi de la transformation de LapLaceE ne figurent pas toutes 
dans le probléme donné [3], on utilise une transformation de FouRIER. 

_ Les développements rappelés dans le premier paragraphe rendent ces distinc- 
tions inutiles et permettent d’opérer wniformément. 

On transforme d’abord le probléme proposé par la transformation de LAPLACE; 


on écrit ensuite immédiatement, sous forme dune-série de fonctions propres, la: 


solution du nouveau probléme; la solution cherchée s’en déduit facilement. 

Outre la simplicité et la rapidité, les principaux avantages de cette méthode 
sont les suivants: 

a) Elle convient parfaitement a |’intégration des équations dont le second 
membre est une fonction génévalement continue — de telles équations se présentent 
fréquemment dans les applications — contrairement aux méthodes qui exigent le 
calcul d’une intégrale générale. 

b) Aucun recours a la théorie des résidus d’ot absence des difficultés éven- 
tuelles signalées plus haut; on reste dans le domaine des variables réelles. 

c) Il n’est pas nécessaire de réduire préalablement les problémes comple- 
tement non homogénes 4 des problémes semi-homogénes. 


3. Exemples 


Hes; deux exemples que nous avons choisis présentent chacun une particu- 
larité: dans le premier, l’équation transformée contient une fonction non donnée 
dans le probléme proposé; les fonctions propres eee interviennent dans le second 
exemple ne sont pas orthogonales. 

On trouvera dans [1] une application a |’équation de la chaleur. 


2) La méthode peut étre étendue au cas de plus de deux variables indépendantes. 
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300 
a) Soit U(x, vy) un potentiel ou une température indépendante du temps 
(équilibre thermique) satisfaisant a l’équation et aux conditions suivantes: 


o?U oU 
rere 


Pet!) ay G(7) (0<*</), 
|U(4,y)| <M (0S% Shy 209), 


Jim [or Ula, 9) + 2 Uz(% ¥)] = A, | 
(y > 0) 


jim [Bs U(#, 9) + Pi Ul = B | 


ot A, B et M sont des constantes; (a? + «3) (83 + 63) + 0, 4% <9, B; By 2 0 


et «1, 8; non tous deux nuls, par hypothése. 
L’application de la transformation de Laprace donne, si l’on pose 


2{ U(x, y)} = ula, s) et H(x) = lim ,U,(x, 9), 
u" + s?u=s5 G(*) + H(*) + que 


a, u(0, s) + a, w’(0, s) = ra 


B, u(t, s) + Biv’, s) = 2s 
Ona (1) 


= u {x 1 
POR SKS \N, (4 ~ 9, B) — 5 G,— Hy — 


X (a, SIN M, ¥ — & Wz COS, %X); 


@,, N;, et 0;, voir page 298; 
1 


ik : 
= x, [TO (a, sin, € — a, w, cosm, &)d& (T=F,G, H). 
0 


Le probléme transformé ne peut avoir une solution pour s = w, (k=1, 2,...) que st 


{ Ox y 1 F; 
— (A — 9, B) ——sG,— H, — — = 
\N, ( Cn ) 5 SG k s ie 0. 
Cette condition détermine H,, (non donné): 
A—g,B ia} 
Te Qk =. Sk 
k N; Ox, G;, On ? 


de sorte que 


_*< i {( @x Wx — S 
U Di eae toe (4 — 9 B) ~ F) SBS + (wy = 8) Gy} 


X (a% SIN MW, ¥ — a Wz, COS Wy *) 


ae che ee pelts 
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pies (x, = 3i((22 (A 2 ore ORY a \ 
bay Wn, — @% B) — :| Oe a te bry, OOK 
1 Oy J 


X (a1 SIN Wy, ¥ — Oy WR COS WR x) . 


Cas particulier [voir aussi (3)] 


a&=fP,;=1, e=—f,=0; A=0; F(x) =0; G(x) = G (constante). 


Cc 
kn : Pe ra all OG alae fen Ve 
ee =D) (=P B : esata eae ates (aU gahelioban 5 eae 
1 G +) : FL U 
fica 1 
Comme 
20 wy (—1)F-1 
¥ = a — sink ] eas 
a sin be x 
a kes oy t= Leet 
arctg > Fe Ah a® sink + Fa. (eal SSM 
Loa cos—- x 1 
ae 
2asin— *¥ co 
l 1 — (—1)* : I 
arctg — a*sink—% (|a¢| S51), 
1— qa? pa k 1 
on a aussi 
CL) ie a el SHG TN AS 
nee 2B e sin i % 26 2a sin x 
U(%, y) i = arctg als arctg Se 


b) Vibrations longitudinales d’une barre de longueur / et de section constante 4. 
L’extrémité x = 0 est fixe; une masse m est concentrée a l’autre extrémité. Ini- 
tialement, la barre est au repos sous l’action d’une force de traction F = AT et 
a Vinstant ¢ = 0, on supprime brusquement F [4]. 

Le déplacement longitudinal U(x, ¢) d’une section transversale de la barre, 
pendant la vibration, satisfait a 


eu, eu rb 

Me faba (2 ar 
ik 

U(x, 0) = ae U;,(#, 0) = 0, 


U(0,t)=0, mU,,(l,t)=—EA U,(i, 2). 


L’application de la transformation de Laprace donne, si & { GME A AES) 


. GO Ss Tx 
IE EA IM OT 
mTl 
u(0,s)=0, ms*u(l,s)+ EA w'(l, A nae 


302 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


Le probléme associé [1] est 
u"+o2u=0, ul0,s)=0, ms* ul, s)+EAu'(l, s) = 0; 


les conditions aux limites contiennent le pavameétre s. 
On liera w et s par l’égalité wm? = — s?/a? de sorte que l’équation caractéristique 
(2) devient 
a2 


cotga@/l= a oO. (4) 


Les fonctions propres sin a, « ne sont pas orthogonales au sens habituel du terme 
mais 
i) 


[sin ow, €sinw,;, dE + 


0 


sinw,lsino,l=0 (i+7). 


La théorie établie dans le cas des fonctions orthogonales [1] reste valable a 
condition de remplacer les intégrales donnant a, et N, par les imtégrales de 
STIELTJES: 

1 


1 : m a : 
a) = a0) fo sin w, € d& + as y(l) sin w;, 1} , 
LO 
1 
i 2 
ns) oh sin?w, €d + 4 sin?a, 1. 
0 
On a (1) 
l 
Ate joer og ae 
a) . vO) | m1 oy Ox + ge | $ sino, $8 sin wz, *¥ 
k 0 
ou 
l sin2 w,/ WO? l sin2 w;, 1 
AAS) nes ke 2 as k 
k 2 Lite at ed oe ae 


en vertu de l’équation caractéristique (4), et 


_ COS Mz d 
ee = m a* we * 
Par suite 
4T sin w, / 
(OER, ve 5 Se ug : 
(+, t) E 4! oy (2 o,1 + sin2 om )) COSA WM, t SIN MW, % « 
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Summary 


Laplace transformation does not seem to be well adapted to the solving of 
certain boundary problems in partial differential equations for two reasons: 
(1) the transformed equation involves a function which is not given in the problem, 
(2) the inversion leads to rather complicated contour integrals (Ale 

We show that even in these cases Laplace transformation is suitable: thanks 
to the results obtained in a precedent work [1] and partially recalled in the first 
paragraph of this paper, we find the solution very easily and quickly; we remain 
in the real domain. 


(Regu: le 26 novembre 1956.) 


Einige Bemerkungen zu M. Hanin und M. Reiner: 
‘On Isotropic Tensor-Functions and the Measure of Deformation’) 


Von HANSWALTER GIESEKUS, Leverkusen, Deutschland?) 


Zu 5: Das von HANIN und REINER entwickelte und in Abschnitt 5 zur Ablei- 
tung einer Tensorformel fiir das Hencky-Mass benutzte Verfahren ist zwar von 
grosser Allgemeinheit, doch kann man in dem vorliegenden Fall den Beweis der 
Formel (6. 1) wesentlich einfacher auf direktem Wege fiihren. Man geht dazu 
von der fiir symmetrische Tensoren mit positiven Hauptwerten eindeutig defi- 
nierten Formel 


a 
In A= /a-*-da 
1 


aus, wobei das Integral als Parameterintegral aufzufassen ist?), und formt hierin 
den Integranden mit Hilfe der Cayley-Hamilton-Gleichung um zu 


1 I II 
ee ye ee gl, | eae 

- ie 4 iA ee 

so dass sich 
| PI oe 
1 2; oe. a / pe | 
iM Arie Jan keh Mies fa 
1 1 1 


ergibt. Setzt man nun speziell 


A 
a=1-2¢de, 
A 


wobei fiir die Invarianten I, II und III von 2 e die folgenden Beziehungen gelten: 


ES as aia al ake lt 


tort Mil a no tol + All, T,=1—ol 4+ el — 6? 11, 


1) Siehe ZAMP 7, 377=93 (1956). 

2) Farbenfabriken Bayer AG. err 

3) Der Kiirze halber verwenden wir hier die indizesfreie symbolische Schreibweise, bezeichnen 
Tensoren durch halbfette Buchstaben und das innere Produkt zweier Tensoren durch einen Punkt. 
Im iibrigen schliessen wir uns jedoch ganz der Bezeichnungsweise von HANIN und REINER an. 
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und integriert iiber den skalaren Parameter o von 0 bis 1, so folgt wegen 


H A 
li EG Ea RAC ee oge'— in (nee) 
unmittelbar 


u 2 1 A 
pap. ted aes) sg 2é-(1- 202) (3-61) , 
carey Han eating QE A) Hees LI GN lool dee Le wn 
h le 
ee ACOs) = 2:0 1 oe oa 
a 1 oho? Tho Ti 
0 


0 


A : ; 
Nach Ausmultiplizieren der Klammern, Umwandlung von (2 Se mit Hilfe der 


Cayley-Hamilton-Gleichung und Ordnen nach Potenzen von 2é erhalt man 
daraus dann leicht die Gleichung (6.1) von HANIN und REINER. 

Man kann diese Formeln fiir die praktische Anwendung noch brauchbarer 
machen, wenn man die Integrale in allgemeiner Form auswertet. Das erfordert 
eine Partialbruchzerlegung und setzt somit die Kenntnis der Wurzeln des allen 
drei Integralen gemeinsamen Nennerpolynoms voraus. Das Problem der Be- 
rechnung dieser Wurzeln ist aber der Bestimmung der Hauptwerte ¢,, ¢, und e; 


A F : 
des Tensors e 4quivalent, denn es gilt: 


iL 1 a 
= oso? Loe Lie = 
ol+o foe 1 Ill ( ea (o =A (o Farle 


Damit erhalt man nach dem bekannten Verfahren im Falle, dass alle Haupt- 
werte verschieden gross sind, fiir die drei Koeffizienten des Tensor-Polynoms: 


€5 € €,@ 
Fi= oes iii ty <=>. Ne ae in(t 2 
E (€1 — &g) (@; — @3) ( Js (€2 — €;) (€2 — @3) ( ») 
, €1 ls 


7 (€3 — €1) (€3 — @g) 


B=3 are Soe ay ml = 2a) + i zy int — 24) 
ae fae ay m2 es), 

a : | a a ae cara In(1 — 2 e,) + c= ae wey. In(1 — 2.) 
eee oe pgeln i2 e.)]- 


Sind zwei Hauptwerte gleich, zum Beispiel e, = es, so ergibt sich statt dessen: 


es 2, — e,) 2 
ky = — a (=e é, (2 ¢ V in(d & & 
. (é. — &)* (¢1 — @)* A ee (€, — @) (1 — 2 e,) ’ 
é 

(een a DS ee 2 in(1 Sea 

PCE amir ree 2 <I ep 

1 1 

= —; In(l — 2e,) 4 inl — & L 

’ 4 (e — ey? | 1) Taleo eee ot) ie ie Bae ae 
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Pur den, physikalisch gesehen, allerdings trivialen Fall, dass alle drei Hauptwerte 
gleich sind, findet man, wie nur aus Griinden der Vollstandigkeit mitgeteilt sei: 
6e?7 —2e L—4e 1 
Pe = = In (f'— 26) fe Re ee 
% i) (1 — 2e,)? ee a a eee ea 
Man erkennt aus dieser Entwicklung deutlich, worin der Vorteil des Hanin- 
Reinerschen Verfahrens vor dem klassischen Verfahren zur Berechnung von 


, ‘ : A 
Tensorfunktionen besteht: Zwar muss man nach wie vor die Hauptwerte von e 
bestimmen, doch erspart man sich die im allgemeinen ziemlich miihsame Aus- 


s é A : : 
rechnung der Transformationsmatrix, welche e auf Diagonalform bringt. 


Zu 7: Der Deformationsgeschwindigkeits-Tensor /,, = (v,,, + Us,,)/2 stellt ein 
stvenges, und zwar im wesentlichen das einzige transformationsinvariante Mass fiir 
die Deformationsgeschwindigkeit dar und ist keineswegs nur durch « premature 
linearization» zu gewinnen. Dieser Tensor beschreibt namlich genau die relative 
zeitliche Anderung des infinitesimalen Abstandes zweier kérperfester Punkte. Im 
Gegensatz dazu ist die zeitliche Ableitung des Henckyschen Deformationstensors 
im allgemeinen kein Mass fiir die Deformationsgeschwindigkeit. Nur im Falle 
einer rotationsfreien Str6mung werden — wenn man noch in bestimmter Weise 
tiber das Koordinatensystem verfiigt — beide Tensoren identisch. Letzteres gilt 
zum Beispiel fiir das von HaniIn und REINER auf Seite 391 angefiihrte Beispiel 
eines eindimensionalen Str6mungsvorganges, wie man anhand der dort gege- 
benen Formeln leicht verifizieren kann. 

Wegen einer allgemeinen Begriindung dieser Behauptungen muss auf die 
Literatur verwiesen werden (vgl. zum Beispiel OLDROyD [1]*), RrvLin und ErRIck- 
SEN [2], sowie TRUESDELL [3]). Doch soll an dieser Stelle der Sachverhalt wenig- 
stens an einem einfachen Beispiel anschaulich gemacht werden. Wir wahlen daftr 
das in obiger Arbeit auf den Seiten 389-390 behandelte Beispiel der einfachen 
Schubdeformation («simple shear») und setzen darin S = gt ein. Dann liegt der 
Fall einer stationaren Scherstr6mung vor. Der Deformationsgeschwindigkeits- 
Tensor lautet dafiir: 


Oo Oo of 
q || | 

hee 7 Wo Oo “inh, 
10 2 20 


wahrend sich fiir die zeitliche Ableitung des Hencky-Masses gemass Gleichung 
(6.13) ein Ausdruck von der Form 


Sar Oar Ap ore lh eens De a 
oe ee ae o if/+4 fio 1 ojft+--- 
Lirara dled dG 0 9 -1|| 


ergibt. Hier stimmt der erste Term mit /,, tiberein, alle folgenden jedoch enthalten 
die Zeit ¢ explizit. Da aber fiir einen stationéren Vorgang ¢ nur bis auf eine will- 
kiirliche Konstante bestimmt ist, so kann dé,,/df unmdglich zur Beschreibung 


eines solchen Vorganges dienen. 
Man kann den diesbeziiglichen Uberlegungen von HANIN und REINER nur so 


einen Sinn abgewinnen, dass man in obigem Ausdruck anstelle von ¢ eine Kon- 


stante mit der Dimension einer Zeit t, das heisst eine Art Relaxationszeit ein- 


4) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 306. 
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fiihrt. Anschaulich liegt dem dann folgende Hypothese zugrunde: Das betrach- 
tete Material soll bei einer stationaren Strémung (bei der Deformations- und 
Riickbildungsvorgange im Gleichgewicht sind) das gleiche Spannungsverhalten 
zeigen wie im Augenblick einer momentan ausgetibten Deformation von der 
Grésse S =qt. Dadurch wird indessen nur eine recht spezielle Klasse von 
Fliissigkeiten erfasst; sie scheint wenigstens in erster Naherung mit der von 
Fromm [4] und prWirt [5] postulierten verallgemeinerten Maxwell-Flissigkeit 
iibereinzustimmen. Solche Fliissigkeiten haben in der Tat keine Querviskositat 
(«cross viscosity »), sind aber bekanntlich keineswegs unelastisch. Auch zeigen sie 
nicht den beschriebenen Zentripetal-Pump-Effekt, sondern vielmehr einen Zentri- 
fugal-Effekt. Dies folgt aus [6], wo gezeigt wird, dass fiir das Auftreten des Zentri- 
petal-Pump-Effekts (dort Torsions-Quereffekt genannt) eine Ungleichung der 
Form 3 Ne - t, > 0 erfiillt sein muss, wobei Ny den Grenzwert der Querviskositat 


und ie denjenigen der ersten « Drehviskositats-Funktion » bedeuten. Diese Gréssen 
sind fiir cine Fromm-deWitt-Flissigkeit mit der « Viskositat» w, aber durch 


ee = 0, ¢, = mT gegeben, so dass obiger Ausdruck negativ wird. 
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Summary 


A simplified proof is given of the Hanin-Reiner formula for the Hencky 
measure of deformation, and the integrals of this formula are calculated. In the 
second part the author shows that, for the particular case of simple shear flow, 
the rate of the Hencky measure is not suitable for the description of the rate of 
deformation. 


(Eingegangen: 8. Januar 1957.) 


Comments on the Preceding Paper 


We are grateful to Dr. GiEsEkus for his interesting remarks. Regarding § 5, 


we appreciate his simple derivation of our equation (6. 1). However, with regard - 


to his further development, if the principal values e,, ¢,, e, of the tensor are 
known (but they are in general not known), all other functions of this tensor 
are also known, and therefore all measures can be defined. The detailed evaluation 
of the integrals, therefore, does not help. 


With regard to § 7, we still regard the classical f, ,, as one possible flow tensor, 


but not the only one. Actually if f,, is an invariant measure, any properly con- 
structed function of /,, will also be a measure. M. Hann and M. REINER 


(Received: April 15, 1957.) 
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Note on a Direct Method of Solving Problems of Elastic Plates with 
Circular Boundaries Having Prescribed Displacement 


By BIBHUTIBHUSAN SEN, Calcutta, India?) 


Introduction 


In a previous paper |2]*), a direct method of solving the first boundary value 
problem of elasticity connected with plates having a straight or circular boundary 
was discussed by the author. The object of this note is to deduce a direct method 
of solving the corresponding problems of circular holes in an elastic plate of iso- 
tropic material with displacements pre-scribed on the edge. To illustrate the method 
the problem of deformation in an infinite isotropic plate due to a nucleus in the 
form of a centre of dilatation at a finite distance from the hole is solved, it being 
assumed that the edge of the hole is rigidly fixed. 


1. Method of Solution 


We take w and v to be the components of displacement along the rectangular 
axes OX, OY, and suppose that they are independent of the coordinate Z. Assum- 
ing further that the Z-component of displacement is zero, we get in the absence 
of body forces the equations of equilibrium 


" 1 ae, 
Eo a ag EE. re ee ns) 
where ‘ ; if a 
uU ) eta. Sal gore 
Og Toy es te Ore Oe. 


and o = Porsson’s ratio (which takes a modified form in the case of generalized 
plane stress). 
A being a harmonic function, we can put 


A= ReG@z) , (1.2) 


in which Re denotes the real part, z= 7 + iy and @(z) is an analytic function 
of z to be determined. 

If the boundary of the circle be v = a, it can be easily found that the solutions 
of the equations (1.1) can be expressed as 


“= aii Re| (a? — r*) 2) + Re[g,(z)] , ae 
. @D 3 
ome eer Re | (a? — r*) a + Re[g,(z)] . 


In the above expressions g,(z) and g,(z) are analytic functions of z such that 


[Regi(2)png = [ly_a> (Regellpna = yaa: (1.4) 


1) Jadavpur University. 
2) Numbers in brackets refer to References, page 309. 


PR. 
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From the given values of w and v on the circular boundary, it is possible in 
general to obtain the functions g,(z) and g,(z). From (1.2) we have 


Ou ov 


— bis oe (1. 5) 


Re[®(z)] = 4d = 
Calculating the values of du/dx and dv/dy from (1.3) we find that the above 
result is satisfied if we take 
1— 2 
@(z) = 2 = 20) 


a a a eae) veel (1.6) 


primes indicating that the functions have been differentiated with respect to z. 
Since g,(z) and g,(z) are supposed to be known, the value of @(z) is given by 
(1.6) and hence from (1.3), w and v are completely determined. 


2. Example 


Suppose in an infinite plate, there is a centre of dilatation in the form of a 
nucleus of thermo-elastic strain at a point A(c, 0), c being greater than the radius 
a of the hole. The boundary of the hole is supposed to be rigidly fixed. Compo- 
nents of displacement due to such a nucleus are [1] 


0 0 
Oe = Mes We USSTal sy Se Ge (log7;) , (2.1) 


where # is a known constant and 
vt = (x — cj? 4+ y?. (2.2) 


Let B be the inverse point of A with respect to the circle, and R,, R,, the distances 
of a point P on the circle from A and B respectively. Then 


6 Cah 6 k Zz 
le ee eee) 
[M],4 R? te pe R3 r=a c es z— a?/¢}\r=a’ 


: at . (2.3) 
hae * Leh. Leh Ge | 


If the circular boundary of the hole be rigidly fixed, we are to find us, v, from 
(1.3) such that i : 


(Milne = — [Milan = = [Re(=—ar)] 


| 


c Z— a*/c 
R 42 Sa 
[Melrna = —[ilag = = C [Re(, = a®/c ipa 
Thus we can put in (1.3) 

_&k ae 1 L ke 2 1 | 
ad=— [1+ Emap], ae--A [1+ 2-45]. eo 
Hence from (1.6) ae eee 

POs c2(3—46)  (¢—a%Jc)?° Cs 


Substituting the values of @(z), g,(z) and g.(z) from (2.6) and (2.5) into (1.3) 
we get uw, and v,. Values of displacement will then be given by 


U=Uy+U,, VH= V4 0,5. (aw) 


a o-e 
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A circular hole with fixed edge in an infinite plate. 
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Zusammenfassung 


Zweck dieser Mitteilung ist die Herleitung einer direkten Methode fiir die 
Lésung des zweiten Randwertproblems in einer elastischen Scheibe mit kreis- 
formiger Begrenzung. Zur Veranschaulichung der Methode wird das Problem der 
unendlichen Scheibe mit einer thermoelastischen Verzerrungsquelle in endlichem 
Abstand von einem kreisfo6rmigen Loch mit fixiertem Rand gelést. 


(Received: February 4, 1957.) 


Elementare Falle des Dirichletschen Problems 
fiir elliptische Gebiete der Ebene 


Von VACLAV VopicKA, Pilsen, CSR 


Es ist nicht allgemein bekannt, dass und wie sich manche interessante Pro- 
bleme der Laplaceschen Gleichung auf ganz elementarem Wege lésen lassen. Das 
zugehorige Verfahren soll hier fiir den Fall spezieller Bedingungen auf der 
Peripherie eines elliptischen Grundgebietes durchgefiihrt werden. 


1. Problemstellung. Es sei f(x, vy) ein Polynom vom beliebigen Grade n = 9. 
Man fragt nach der Lésung u(x, vy) der Potentialaufgabe: 


_ OF O°u 
Ox? oy? ’ 
u(x, y) = f(#, y) auf der Ellipse 2 47 == a? 2 "a2 be = 0: 


Au 


= 0 im Gebiete b? 72+ a? y?—a?b?<0, 


(1) 


2. Lésung. Die Grundlage zur elementaren Behandlung unseres Problems 
bilden zwei bekannte Satze der Algebra, die wir ohne Beweis in einer fiir fol- 
gende Zwecke passenden Form aussprechen wollen. 
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Satz 1. Ein Polynom (x, y) n-ten Grades hat dann und nur dann tiberall auf 
dey Ellipse b? x® + a® y? — a® b? = 0 den Wert Null, wenn 


p(¥, Y) = (0% #2 + a® y? — a* b*) w(x, 9); (2) 
hievbei ist w(x, y) beliebiges Polynom vom Grade n— 2, falls n22 und 
p(s, vy) =0 fir n S11. 

Satz 2. Weiss man, dass ein Polynom (x, vy) n-ten Grades (n = 0) mit der 


ganzen vationalen Funktion b® x® + a? y? — a® b® mindestens 2 + 1 verschiedene 
Nullstellen gemeinsam hat, so ist p(x, y) = 0 uiberall auf der Ellipse 


6b? #2? + a* y*— a®b?=O0, 
und nach Satz 1 ldsst sich dann (x, vy) in der Form (2) ausdricken. 


Mit einer Nullstelle irgendeiner Funktion F(x, y) meinen wir natiirlich jedes 
Zahlenpaar (¥, y) mit der Eigenschaft F(%, y) = 0. 

Nunmehr kehren wir wieder zum Problem (1) zuriick und setzen die Lésung 
u(x, y) mit noch unbekannten 2” + 1 Koeffizienten yo, yy, yy (vy = 1, 2, 3,... , ”) 
in der Form 4 


ul, 9) = Yo +S ty Pl%, ”) + ro golx, WI (3) 


ah 
voraus. Hierin sind 


‘ 


le | 
pul 9) = S)(—1)"(0),) 7 w(*, 9) = (-1"(, ae 1) onto (3.11 


2=0 2=0 


wig 


die durch Soret ; 

(w+ ty)" = polx, vy) +t Q(*, ¥) (3.2) 
definierten harmonischen Grundpolynome, und [vy] bedeutet ganz allgemein den 
ganzen Teil der Zahl 7. 

Die Beiwerte y aus (3) bestimmen sich nun einfach aus dem linearen Glei- 
chungssystem 


n 
U(%5, Vs) = Yo +2 'lv Pil% 5 Vs) + 1 Ge V— Kee d,) (8S Woe See 
y=] 


worin (*,, ¥,); Ss = 1, 2,..., 2+ 1 beliebige untereinander verschiedene Punkte 
der Ellipse b? x? + a y? — a?b? = 0 sind. In allen praktisch wichtigen Fallen 
sind die (¥%,, y,) immer so wahlbar, dass die Beziehungen (4) und damit auch alle 
folgenden Rechnungen verhaltnismdssig einfach ausfallen. 

Es ist leicht einzusehen, dass die Gleichungen (4) lésbar sind, und das Ergebnis 
der Einsetzung der daraus berechneten Werte y in (3) schreibt sich dann in Form 
der Determinantengleichung 


Is Pil*, y), q(%, y), P2(*, y), qo(*, sai boty. Pr l*, y), In(*, y), u(x, y) 

1, Pi(Pr), qi(P,), po(P,), q2(F), =e Palla), Yn(Pi); (Pr) 

1, p,(Ps), qi(Po), P2(Ps), q2(Po), «Seah as InP 2), T(P) 6 
? Pi(Pon); Q1(P2n)s Po(Pen), Jo(Pon), ’ Puli ou) Gatos {(Poe) 

tr Palbenei)) Witenes), Pe(Pon +1): do(Pan4a)s «==> Pa(Pan+1)! Gal Pina)’ Ponsa) 


ae Ye 
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Natiirlich ist hierbei 


p(P.) = po(*., Vs) , qv(P,) = qv Xs, Vs) ’ Ts) = I(%s, V5) . 


Dass die durch (5) dargestellte Funktion w(, y) wirklich die gesuchte Lésung 
angibt, sieht man leicht folgendermassen cin: Erstens ist w(*, y) nach (3) ein in der 
ganzen Ebene harmonisches Polynom n-ten Grades. Weiter haben nach (4) die 
beiden Polynome u(x, y) — f(x, y), b2 x2 + a® y2— a2 b? alle 2n+1 Punkte 
(¥,, Vs) als gemeinsame Nullstellen. Da u(x, y) — f(x, vy) vom héchstens n-ten 
Grade ist, so hat man nach Satz 2 wirklich u(x, y) — f(x, y) = 0 langs der ganzen 
Ellipse 62 #2 + a? y2 — @2b%= 0. 

Nach demselben Satz 2 lasst sich unsere durch (5) definierte Losung auch in 
der Form 

u(x, y) = f(x, vy) + (b? x8 + a? y? — a? b?) w(x, y) (6) 


schreiben. Dabei ist p(x, v) ein passendes Polynom vom héchstens (n -- 2)-ten 
Grade, oder y(¥, vy) = 0. Die genaue Bestimmung von y(x, y) ist in konkreten 
Fallen nicht schwer, wie sich noch weiter zeigen wird. 


3. Tabelle der Polynome (3.1). Die ersten 10 Paare von Polynomen (3.1) sind: 


pPy=* Gis: 

y= **— y* z= 2xy 

Py —3 xy" eS = yg" 

py= x*— 6x2 y?2 4 yt d= +e y—4nr y3 

p= #>— 1045 2 +5 x y* Gg = 5x4 y— 10 x? yF + y5 


Pee #8 — 15 x4 yy? +15 e72 yt — y® Ig= Ori y— 2048 y8 + 6% yd 
py = e7— 21 xi y2+ 35 x8 yt— Trey gg= 7x8 y — 35 wt y3 + 21 x2 yS— y? | (7) 
ps = 8 — 28 x8 y2 + 70 x4 y4*— 28 x? yS + y8 
: Ig= 8x? y— 56x > y3 + 56048 y>— Bry? 
Dp, = 2° — 36 x7 y* + 126 x y* — 84 x8 48 + 9% y8 
dg = 9 x8 y — 84 x8 y® +7126 x4 y® — 36 2 y7 + 9? 
Pip = #19 — 45 x8 y? + 210 x8 y* — 210 x4 y8 + 45 4? y8 — y 
Gp = 10 %® y — 120 #7 y? + 252 #5 y? = 120 x8 yy? + 10 4% 9? 
4. Beispiele. Es soll jetzt obige Theorie auf einigen Beispielen erlautert werden. 
4. 1. Ist f(x, y) selbst ein harmonisches Polynom, so folgt aus (6) fiir p(¥, y) 
die Bedingung Get et ott ave Te. 


welche im ganzen Gebiet b? x? + a® y? — a? b2 < 0 erfiillt sein soll. Dies tritt dann 
und nur dann ein, wenn das Polynom in der eckigen Klammer identisch ver- 
schwindet, das heisst p(x, y) = 0. 

Im Sinne von (6) kommt man so zum folgenden (hier freilich fast trivialen) 
Ergebnis: Ist f(x, y) beliebiges harmonisches Polynom, so lautet die Lésung 
u(x, vy) des Problems (1) ganz einfach 


u(x, vy) = f(%, 9) - 
Insbesondere gehoéren hierher die haufigen Falle . 
t(%, ¥) = foo: f(¥, ¥) = foo + fro * + for ¥ 


mit beliebigen konstanten f; ,. 
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Im weiteren kann man sich also auf nichtharmonische Polynome /(*%, y) be- 
schranken. 
4, 2. Es sei 
Hai v= foo F higt a fo VE fee a foo 75 


fin gegebene Festwerte. 
Mit den 5 Punkten 


(4 @,0), (0, £9), (eae = 


der Randellipse b2 x? + a? y2 — a? b? = 0 lautet die Gleichung (5) [siehe auch (7)]: 


1, eh 05 VES, EY, UD) 

LS ae Oras, 0, foot 4fhio + 2? foo 

Lae 0, “<a, 9, foo — 4 fio + 4? foo | ee fa 
Lire (Oimgs Oy made 0, foo + © for + 5? fos 
ES RUSS alae 0, foo — 8 for + 2? fos 

iL, a, a, 0, a, foook & Vig toner (foo + far + foe) 


Mit einiger Rechenarbeit ergibt sich daraus die nicht allgemein bekannte 
Losung a pt 
a 


a? + db? (feo + foo) + fro ¥ + for Y 


u(x, ¥) = foo + 

(8) 

a* fog — 5? foo (a2 
a? + b? ! 


ale y)y+firy. 


Fur den Spezialfall b = a erhalt man die Lésungsform 


u(x, Y) = foo + a (foo + foo) + fro ¥ + for ¥ + E (foo — fos) (47 — 7) + fave (84) 


des Dirichletschen Problems fiir den Kreis mit dem Radius a. 


Bemerkung. Das Ergebnis (8) folgt in ein wenig anderer Gestalt 
u(%, ¥) = foo + fio ¥ + for ¥ + foo ¥* + fir ¥ ¥ + Soo Y?” 
9 
-; ine : fe (0? 42 + a? y? — a2 62) BL 
viel schneller aus (6). 
4, 3. Fir den Fall 
f(#, ¥) = foo + fio * + for Y¥ + foo ¥® + fr ¥ Y + fos V? 
+ foo #8 + for 47 VY + fie ¥ Y® + fos Y® - 
fiihrt der Ansatz (6), das heisst 
w(x, ¥) = T(¥, y) + (D? 4? + a? y® — a? b*) (yy Cy % + Opi Y) 
zur Bedingung 
foo + fos + (3 fso + fae) ¥ + (fer + 3 fos) ¥ + 2 (Cy 02 % + Cox 2 Y) 
+ (a? + b?) (Coo + Cig ¥ + Cp, Y) = 0, 
die im ganzen Gebiete b? x2 + a? y? — a2 b?< 0 erfiillt sein soll. 
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Dies ist der Fall, wenn und nur wenn 
Coo (@* + B*) + fon + foo = 0, Crp (47. +3 B72) + fy, + 3 fo = 0, 
Coy (3. a* + 5) + 3 fog fey = 0. 


Mit den daraus berechneten Werten von c,, ergibt sich dann die Lésung 


u(x, Y) = foo + fio * + for ¥ + fay ¥*° + far VY + foo v? 
+ feo ¥* + for *7Y + fie ¥ ¥? + fos ¥® | (10) 
feo + fos fie > 3 tao 3 fos + fe 


Technisch sind unsere Ausfiihrungen iiberall dort von Interesse (Elektro- 
statik, Warmeleitung usw.), wo man mit der Laplaceschen Gleichung zu tun hat. 

Es ist zu bemerken, dass sich auch andere Fragen der mathematischen Physik, 
insbesondere unter gewissen Umstanden das Poissonsche Problem 


Au = f(*, y), w= 0 am Rande des Grundgebietes 


nach der vorangehenden Methode lésen lassen. 


Summary 


The Dirichlet problem is solved for an elliptic domain if the boundary values 
are given by a polynomial defined in the interior of the ellipse. The solution is 
a finite expansion in harmonic polynomials. 


(Eingegangen: 17. Dezember 1956.) 
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Friihjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 4. bis 5. Mai 1957 in Brunnen 


Berichte iiber die Sitzungen fiir angewandte Physik und Mathematik 


Ein Widerstandsbrett fiir rotationssymmetrische Potentialprobleme. 
Von E. BAUMGARTNER, R. Gatti und P. HuBER, Basel?). 

Seit der praktischen Realisierung von Widerstandsnetzwerken durch DE 
PackH?) und RepsHaw®) wurde diese Art Analogierechengerat auch auf kompli- 
ziertere Differentialgleichungen mit Erfolg angewendet*). Diese Methode bietet 
gegeniiber dem elektrolytischen Trog hauptsachlich den Vorteil grosserer Genauig- 
keit und einfacherer Handhabung. Ein scheinbarer Nachteil tritt aber bei 
krummen Oberflachen auf, weil ohne jeweilige sso der Anordnung nur 
diskrete exetenpensic zur Verfiigung stehen. 


1) Physikalisches Institut der Universitat. 
) D. C. pe Pacxu, Rev. sci. Instrum. 18, 798 (1947). 
3) S.C. Repsuaw, Proc. Instn. mech. Engrs., Lond. 159, 25 (1948). 
4) C, Lrepmann, J. int. Calc. analog. 1955, 346. 


bo 
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In Figur 1 ist neben dem Prinzipschema eines solchen Widerstandsnetzwerkes 
dessen praktische Realisierung dargestellt. Bei einem bestimmten Problem wer- 
den entsprechend den Randbedingungen Potentiale auf das Brett aufgedriickt, 
dann stellt sich in den anderen Punkten eine angenaherte Losung der Laplace- 
Gleichung ein. Alle Widerstande des Netzwerkes wurden mit Konstantandraht 
auf hochohmigen Widersténden als Wickelkérper gewickelt, und zwar auf 
+ 1,5°/)) genau. Eine geeignete Wicklungsart erlaubt die Beniitzung von Wech- 


Figur 1° 
Widerstandsbrett und Prinzipschema des/Netzwerkes. In der radialen v-Richtung sind 11, in der 
axialen z-Richtung 40 Maschen. R = 8k. 


“e 


oe 
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selstrom (50 Hz), ohne dass Schwierigkeiten wegen Blindwiderstanden auftreten. 
Die Lésung*) der allgemeineren Gleichung div (e grad V) = —9 erlautert Figur 2. 


Figur 2 


Berticksichtigung von 
a) Raumladungen durch die Zufiihrung von Strémen. 


oA ¢ +0); 7 = —— oA? =). 


Hierbei ist @ jeweils der Mittelwert von g im eingezeichneten Volumen v. 


b) Dielektrizitatskonstanten, jeweils durch Ersetzen des Widerstandes R; durch R,/é, wobei € der 
Mittelwert von € auf der Oberflache o ist. 


491 Tan 627 744 


° 
626 742 


333 $56 714 833 926 


Figur 3 
Potentialverteilung von zwei konzentrischen Kugeln. Die Pfeile bedeuten Stromzufiihrungen. Die 
oberen Zahlen sind aus gemessenen Werten ermittelt, die unteren sind die theoretischen Sollwerte. 


5) C. Lresmann, Brit. J. appl. Physics, 7, 92 (1956). 
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Randbedingungen, die nicht auf Knotenpunkte fallen, werden wesentlich 
rascher durch Zusatzstroéme realisiert als durch Widerstandsabanderung. Durch 
geniigend hochohmige Potentiometer (auf das Brett aufsteckbar), deren Abgriff 
auf den Berandungsort eingestellt wird, lasst sich leicht der notwendige Strom 
einstellen. In Figur 3 wird dies durch ein Testproblem illustriert. In diesem wurde 
der Maschenfehler durch Ausmessen des Problems in einem halb so grossen 
Ma8stab weitgehend eliminiert. 

Das Brett wird jetzt fiir ionenoptische Probleme im Zusammenhang mit dem 
Bau eines Beschleunigers verwendet. 


Zur Theorie der Fligelgitter. Von L. Meyer, Ziirich"). 

Seit langerer Zeit beschaftigt sich das Institut fiir Aerodynamik der ETH mit 
der Berechnung von Schaufelprofilen fiir axiale Turbomaschinen. Da bei den 
Rechnungen von reibungsfreier Stromung (Potentialstromung) ausgegangen 
wird, werden die Experimente zur Priifung der Theorie mit der Methode der 
elektrolytischen Analogie durchgefiihrt. Der ebene Trog enthalt eine endliche 
Anzahl Profile des ebenen Gitters; als Berandungen werden je zwei Stromlinien 
und eine Potentiallinie der An- und Abstromung materialisiert. Der aerodyna- 
mischen Zirkulation jeder Schaufel muss eine elektrische Periodizitat entsprechen: 
man verwendet daher die konjugierte Analogie 


il 
Dio ae tae gee agit o/h ®,, 
(o = spez. Widerstand, h = Hohe der leitenden Schicht). Der Zirkulation IF ent- 
spricht dann der Strom J durch die Schaufeln; diese, gespeist durch Potentio- 
meter, sind als Elektroden ausgebildet, ebenso die die Stromlinien darstellenden 
Wande, wahrend die senkrechten Abschliisse der letzteren aus Isoliermaterial 


bestehen (Figur 1). Die Messungen bezwecken die Priifung des Umlenkverhaltens 


a= V-aN +aV a= V,- 4,V Vy 
\ 


Figur 1 
Ebenes Schaufelgitter im Elektrolyttrog. 


1) Institut fiir Aerodynamik der ETH. 


eee Ey 


Se 
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(Abstr6mwinkel) des Gitters bei verschiedenen Anstrémwinkeln sowie die Kon- 
trolle der berechneten Geschwindigkeitsverteilungen am Profil. 

Es hat sich gezeigt, dass neben grosser Genanigkeit der geometrischen Anord- 
nung eine sorgfaltige Einstellung der Randbedingungen des elektrischen Feldes 
notwendig ist. Die Periodizitat verlangt gleiche Spannungsdifferenzen von 
Schaufel zu Schaufel, sowie gleiche’ Stréme J durch alle freistehenden Schaufeln. 
Sind f, und f, die Zu- und Abstrémwinkel in bezug auf die Gitterrichtung, so 
bedingt die Analogie von J’und J ferner, dass 


i “ajh (N+ 1) | (ctg B, — ctg B,) | = Jen 


sein sollte (N = Anzahl freistehender Schaufeln, hier zum Beispiel 5; o/h wurde 
wahrend des Versuchs mit einem aus 2 koaxialen Zylindern bestehenden, in?) 
beschriebenen Apparat gemessen). Figur 2 zeigt links Spannung und Strom der 


B, =182° 


\aV 
sy=0972 


| 


eS Oy ay = 5, AV AV; 


| 


LF iis ee = 0015; AV =AV,,:Js = Jeni J1,3,4,5™ Jin - 
Vo ~ M8) 
Figur 2 
Links Spannung und Strom der Schaufeln, rechts Geschwindigkeitsverteilung entlang der mitt- 
leren Schaufel fiir 3 Stadien der Einstellung des elektrischen Feldes. 


2) P. pE Hatter, Anwendung der elektrischen Analogie zur Untersuchung von Schaufelgittern, 
Techn. Rdsch. Sulzer 1947, Nr. 3/4. 
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Schaufeln, rechts die auf Grund der Analogie mit einer Sonde gemessene Ge- 
schwindigkeitsverteilung’) an der mittleren Schaufel fiir drei Stadien der Ein- 
stellung des elektrischen Feldes. Bei den Einstellungen @ und @ wurden die 
aerodynamischen Stromlinien der Randschaufeln mit den bei deren Nase und 
Hinterkante beginnenden geraden Begrenzungsplatten identifiziert. Nach Erfiil- 
lung der Periodizitatsbedingung fiir die Spannung (Einstellung @) ergab sich 
zunachst ungeniigende Periodizitat und Grésse der Schaufelstr6me; an der mitt- 
leren Schaufel wurde etwa 6% zu wenig Strom gemessen, welches Manko sich in 
einem zu kleinen Integral der tiber der abgewickelten Profilkontur aufgetragenen 
Geschwindigkeitsverteilung aussert. Bei Einstellung @® wurden die Stréme durch 
alle freistehenden Schaufeln auf den geforderten Betrag einreguliert; hierbei 
musste auf die strenge Periodizitat der Schaufelspannungen verzichtet werden, 
Bei Einstellung ® wurde nun beriicksichtigt, dass die Schaufelstromlinien wegen 
ihrer Kriimmung in Schaufelnahe auch in grésserer Entfernung von den Schau- 
feln im allgemeinen nicht mit den zu ihnen parallelen Geraden durch Nase und 
Hinterkante der Profile zusammenfallen, sondern gegentiber diesen etwas ver- 
schoben sind. Diese Verschiebung wurde fiir die Stromlinien der Randschaufeln 
durchgefiihrt, wobei von den beiden Moglichkeiten der mechanischen oder elek- 
trischen Verschiebung die letztere einfacher zu realisieren war (siehe Figur 1). 
Wahrend die Theorie fiir das vorliegende Gitter eine vernachlassigbar kleine Ver- 
schiebung A,V der Stromlinien der Abstro6mung ergab, wurde die mit groésserem 
Anstellwinkel zunehmende Verschiebung 4,V jeweils so eingestellt, dass bei 
strenger Periodizitat der Schaufelspannungen die mittlere Schaufel vom richtigen 
Strom J;, durchflossen wurde. Aus Figur 2 geht hervor, dass dann auch die 
Stré6me durch die iibrigen Schaufeln, insbesondere durch die der mittleren 
Schaufel benachbarten, den geforderten Wert bis auf einen kleinen Bruchteil 
erreichten. 


Farbmetrische Kennzeichnung von Farbstoffen. Von E. Ganz, Basel?). 

Das farbmetrische System der internationalen Beleuchtungskommission (IBK) 
hat sich zur objektiven Kennzeichnung von Farbungen, Anstrichen und farbigen 
Lichtern bewahrt; es eignet sich aber nicht zur Charakterisierung von Farbstof- 
fen, da dessen Koordinaten nicht den farberischen Begriffen von Farbstarke und 
Reinheit eines Farbstoffes entsprechen. Die Darstellung der Farborte der Aus- 
farbungen eines Farbstoffes in zunehmender Konzentration ergibt Kurven, die 
vom nahezu weissen, ungefarbten Substrat bei stetig abnehmender Helligkeit bis 
zu einem Sattigungsmaximum fiihrt und bei weiterer Steigerung der Farbstoff- 
konzentration mit abnehmender Sattigung gegen den Schwarzpunkt lauft [1]?). 
Diese Kurven gleicher Farbstoffreinheit stellen mit den Kurven gleicher farbe- 
rischer Starke das zur Farbstoffkennzeichnung geeignete Koordinatensystem 
dar. Wir haben versucht, dieses Farbsystem aus dem System der IBK auf Grund 
folgender einfacher Idealisierungen der physikalischen Eigenschaften von Farb- 
stoff und Ausfarbung zu berechnen: 


1. Der Verlauf der Extinktion e, der farbbedingenden Absorptionsbande wird 
durch eine Gaufsche Verteilungsfunktion mit der Bandenmitte 4, und der Halb- 


3) Die Bezugsgeschwindigkeit V, (= V4) ist hier die zur Gitterrichtung normale Komponente 
der Eintrittsgeschwindigkeit Cj. : 

1) CIBA AG. 

*) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 323. 
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wertsbreite 6 dargestellt: 
ay Ajo 


& = €° ais 

2. Die Remission R der Farbungen als Funktion von Extinktion e, und Farb- 

stoffkonzentration ¢ wird mit einer modifizierten Formel nach KuUBELKA-MuUNK [2]- 

‘PINEO [3] dargestellt, wobei 7 die Oberflachenremission des schwarzgefarbten Sub- 
strates bedeutet: 


B= i — ae) (UR eee ok); Hy ac tp. 
y wurde hier (fiir Textilien) = 0,02 gesetzt [4]. 


3. Die IBK-Farbwerte dieser idealisierten Farbungen wurden mit einem 
IBM-Rechenlocher 604, Modell 4, durch Summation bei einem Wellenlangen- 
inkrement von 3 my fiir 21 iiber das sichtbare Spektrum verteilte Banden (A,) 
bei je 10 Halbwertsbreiten (8) im Bereich von 60 bis 350 my und je 18 iiber einen 
Bereich von 5- 10° gestuften Farbstoffkonzentrationen (c), das heisst 3780mal 
berechnet. Die erhaltenen Farben bedecken den Bereich von Gelb iiber Orange, 
Rot, Purpur, Violett bis Blau. Die Erganzung des fehlenden Bereiches von Blau 
liber Griin zu Gelb, der sich nicht durch eine einfache Absorptionsbande erzeugen 


10 r 
W 


A= 640 mu 


d f 
#-60mp /// 


Wf 


i: ee. EE 
WYP +016 WV) 


’ . . . 
-Idealisierte Farbungen von Farbstoffen verschiedener Halbwertsbreite mit den Bandenmitten 


400 my (Gelb) bzw. 640 mu (Blau). — — — Kurve maximaler Sattigung; —-+—-— Kurve maxi- 
maler Relativ-Sattigung. 
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lasst, wurde sowohl durch analoge Berechnungen mit Doppelbanden (Halb- 
wertsbreite der Liicke gleich Bandenhalbwertsbreite) sowie mit Rechnungen 
komplementarer Farbwerte erreicht. 

Die Schar der erhaltenen Kurven mit der Halbwertsbreite 6 als Parameter 
sind die gesuchten Koordinaten gleicher farberischer Reinheit fiir alle Farbstoffe 
der Bandenmitte 4, (Figur 1). Sie bedecken von Grau mit 6 = oo den ganzen 
moglichen Farbraum bis zu den Schrédingerschen Optimalfarben [5] fiir B > 0. 
Damit wird die Kennzeichnung von Reinheit und Farbton eines Farbstoffes 
durch die Angabe von f und A, méglich. 


4. Zur Kennzeichnung der farberischen, durch Richttypfarbungen definierten 
Farbstarke verbleiben von den 4 Parametern /, f, & und c der idealisierten 
Farbungen noch die beiden letzteren. Es erweist sich als zweckmassig, die Maxi- 
malextinktion ¢, der idealisierten Farbstoffe so zu normieren, dass bei c = 0 


(+514) 
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(544) 
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-3 — = : 
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: Figur 2 
« Chromatic values » idealisierter Farbungen mit umschliessendem Reinstfarbenzug (Optimalfarben — 
grosster Sattigung). . 
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Weiss und ¢ > oo « Schwarz» die Konzentration c= 1 die Richttypstarke K* 
ergibt. Dann ist die reziproke, zur Erreichung der Richttypstarke erforderliche 
Konzentration 1/c* eines realen Farbstoffes gerade die in Einheiten der nor- 
mierten idealisierten Farbstarke e, gemessene Farbstarke e*. 


1 

Bhat he Pe, 5 re ; 

5. Die Richttypstarke ist als empfindungsmetrischer Begriff im Gegensatz zu 
den reizmetrisch fassbaren Begriffen yon farberischer Reinheit und Farbton erst 
nach Transformation der Farbwerte in den empfindungsmetrischen Farbraum 
bestimmbar. Auch friihere Versuche von GopLOoVE [6] und Kocu [7] zur empiri- 
schen Bestimmung der Richttypstarke als Funktion farbmetrischer Gréssen 
stiitzen sich direkt oder indirekt auf empfindungsmetrische Systeme. Als emp- 
findungsmetrische Naherung benutzen wir ADAms «chromatic value» [8] mit Vj 


W-Vx 
3 = 


30 
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Figur 3 

Kurven gleicher Helligkeit Y der Flache grésstméglicher Relativ-Sattigung O,,q- 
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als dritter Koordinate (NIcKERSON [9]), wobei die Munsell-Transformation der 
Farbwerte durch eine Parabel 5, Grades [10] wieder mit dem IBM-Rechenlocher 
und einer IBM-Sortiermaschine durchgefiihrt wurde. Durch Vertauschen von 
Abszisse und Ordinate sowie durch Umkehr der Vorzeichen der beiden Koordi- 
naten ergibt sich gegeniiber ADAms eine Darstellung, deren Orientierung analog 
der Farbtafel der IBK ist (Figur 2). Der mit Farbstoffen einer vorgegebenen 
idealisierten Halbwertsbreite realisierbare Farbbereich lasst sich hier in bezug 
auf die gesdttigtsten Optimalfarben (Reinstfarbenzug) leicht abschatzen. 


6. In diesem angenadhert empfindungsmetrischen Farbraum wurde die Flache 
der Farborte der mit jedem idealisierten Farbstoff erreichbaren grosstmodglichen 
Relativsdttigung o,,,, in bezug auf die entsprechende gesattigte Optimalfarbe 
(Reinstfarbe) graphisch ermittelt (Figur 3). Es ergab sich, dass nur die reinsten 
1/1-Richttypfarbungen [7] in dieser Flache liegen und dass die iibrigen Richt- 
typen um so relativ geringere Helligkeiten besitzen, je triiber die verwendeten 
Farbstoffe sind (Figur 4). Die Arbeiten zur Bestimmung der Richttypflache im 
Farbenraum sind noch im Gang; deren Abweichung gegentiber der Flache grésst- 
moglicher Relativsattigung erlaubt, den koloristisch als Starke bewerteten 
Grau-Anteil abzuschatzen. 


Y Richttyp 


Y6 max, 


ee 


0 Of 02 03 04°) "0s 06 ~ o7 Soaks 
Figur 4 
Relativhelligkeiten der 1/1-Richttypen auf Wolle, nach Kocr [7] in Abhangigkeit der Relativ- 
Sattigung o. 


Die Maschinenrechnungen wurden von der Lochkartenabteilung der CIBA AG, 


Basel (Leitung W. Caspar), durchgefiihrt. Die umfangreichen graphischen Ar- 
beiten besorgte V. FrReETz. 
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Ein universelles Spektralphotometer. Von J. DuRAND und E. Luiscuer, 
Lausanne?). 

Es wurde ein photoelektrisches Spektralphotometer konstruiert fiir einen 
Spektralbereich von 2100 A bis 50000 A. Ein Schema ist in Figur 1 dargestellt. 
Der Monochromator, bestehend aus zwei spharischen Kollimatorspiegeln S, und 
S,, dem Diffraktionsgitter G oder einem Planspiegel S, und einem Prisma P, ist 
in der Aufstellung von CzERNy und TURNER?®). 


Figur 1 


Prinzipschema Spektralphotometer: F,, F, Eintritts- bzw. Austrittsspalt; $1, S, Kollimatoren; 
G Plangitter; S, Planspiegel; P Prisma. 


Nach WrHLm?) erhalt man mit dieser Montage eine Verbreiterung des Bildes 
des Eintrittspaltes um einen Faktor 1,4, bedingt durch den Astigmatismus der 
spharischen Kollimatorspiegel. Das Defilee des Spektrums am Austrittsspalt 
erfolgt durch Rotation des Diffraktionsgitters um eine senkrechte Achse oder 
fiir die Prismenaufstellung durch Rotation des Planspiegels S,,. 


1) Optisches Labor Metrohm, Lausanne. 
2) M. Czerny und A. F. Turner, Z. Phys. 61, 792 (1930). 
3) H. Win, Opt. Acta 2, 50 (1955). 
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Fiir den Wellenlangenbereich 2100-6000 A wird ein Prisma aus homogeni- 
siertem Ultrasil‘), fiir den Bereich 4000-9000 A ein Glasprisma und fiir das nahe 
Infrarot ein NPL-Plangitter®) verwendet. 
Die Breite von Eingangs- und Austrittsspalt kann kontinuierlich zwischen 
10 7 und 2 mm verdndert werden. Die Einstellung erfolgt symmetrisch und simul- 
tan. Der Verstellmechanismus der beiden Spalte ist in Figur 2 dargestellt. 


Figur 2 


Spaltverstellung. 


Als Lichtquellen fiir die obenerwahnten drei Wellenlangenbereiche wurden 
Wasserstoff-, Wolframbandlampe und Nernst-Brenner verwendet. 

Als Strahlungsempfanger kénnen wahlweise Photozellen, Photomultiplier 
oder Thermoelemente verwendet werden. Das Lichtbiindel wird vor dem Ein- 
trittspalt mit einer rotierenden Scheibe moduliert (11 Hz), was eine Wechselspan- 
nungsverstarkung erméglicht und die Empfanger zugleich unempfindlich macht 
fiir nichtmoduliertes Streulicht. 


Uber stabilisierte Netzgerate mit Transistoren. Von E. BALDINGER und 
W. Czaja, Basel?). 

Stabilisierte Netzanschlussgerate fiir niedrige Spannungen und relativ grosse 
Stréme besitzen fiir viele Anwendungen ein grosses Interesse. Zur Konstruktion 
solcher Gerate sind Transistoren wesentlich geeigneter als Rohren. Figur 1 zeigt 
zwei gebrauchliche und von den R6hrenschaltungen her bekannte?) prinzipielle 
Anordnungen. 


Figur 1 
Prinzipielle Stabilisierungsanordnungen. U, = Referenzspannung. 


4) H. Monn, Festschrift Heraeus (Hanau 1951). 

5) L. A. Sayce und G. D. Dew, Proc. roy. Soc. [A] 20, 278 (1951). 

1) Physikalisches Institut der Universitat, Abteilung fiir angewandte Physik. — Eine aus- 
fiihrliche Publikation erscheint demnachst in der ZAMP. 

2) Vergleiche zum Beispiel F. V. Hunt und R.W. Hickman, Rev. sci. Instrum. 10, 6—21 (1939). 
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Es soll am Beispi ‘| é fii Mie rates i 
“ee : be spiel eines ausgefihrten Netzgerdtes gezeigt werden, dass 
urchaus brauchbare Eigenschaften erzielt werden kénnen. In Figur 2 istedas 
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Figur 2 


Stabilisiertes Transistor-Netzgerat. Der Schalter S reduziert bei dauernder geringer Belastung den 
Strom durch die Quertransistoren zur Vermeidung unnétiger Warmeentwicklung. 


Schaltschema eines fiir Laborzwecke hergestellten stabilisierten Transistor-Netz- 
gerates nach dem Prinzip von Figur 1a dargestellt. Die Ausgangsspannung kann 
von 2 V bis 15 V variiert werden, und das Gerat ist bei 2 V bis zu 2,5 A, bei 15 V 
bis zu 0,9 A belastbar. Der Stabilisierungsfaktor 


AU; U4 
ure Uz AU, 


betragt ~ 10%, der Innenwiderstand des Netzgerates R; ~ 0,5 - 1073 2. Der Tem- 
peraturkoeffizient der Ausgangsspannung liegt bei etwa 3 mV/°C. Bei konstanter 
Umgebungstemperatur (Laborbedingungen) bleibt die Ausgangsspannung tiber 
eine Stunde auf etwa + 1 mV konstant. Die Spannungsdifferenz zwischen dem 
kalten und dem warmen Zustand des Gerates betragt etwa 30 mV, und der End- 
zustand wird in 1/2 h nach dem Einschalten praktisch erreicht. 

Als Referenzspannung wird eine stabilisierte Hilfsspannung mit einer Zener- 
Diode verwendet (durch - — - — in Figur 2 angedeutet). In der vorliegenden Aus- 
fiihrung ist die Konstanz der Ausgangsspannung praktisch durch die Eigen- 
schaften dieser Referenzspannung begrenzt. 

Bei hochstabilisierten Netzgeraten besteht ein wesentliches Problem darin, 
unerwiinschte Schwingungen der Regelschaltung zu vermeiden. Zu diesem Zweck 
sind zusatzliche RC-Glieder eingebaut, die den Frequenzgang passend begrenzen. 
Der Innenwiderstand ist bis zu 10 kHz durch die Regulierung gegeben 


5 10-8. 2 Shy 8.1 02.0.4 
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Fiir hohere Frequenzen bestimmen die am Ausgang liegenden Elektrolytkonden- 
satoren den Innenwiderstand, der bis zu 100 kHz 3,5 - 10~* 2 nicht tiberschreitet. 

Fiir die finanzielle Unterstiitzung dieser Untersuchungen sind wir der Stiftung 
Hasler-Werke zu Dank verpflichtet. 


Aufbau grosser Signale mit Elimination von Reflexionen in magneto- 
striktiven Speicherleitungen. Von D. Marper, Ziirich’). 

Speicherleitungen, welche auf der Erzeugung elastischer Wellen in ferroma- 
gnetischen Drahten beruhen, bilden einen wesentlichen Bestandteil der techni- 
schen Realisierung von Vielfachzahlgeraten oder kleinen Rechenmaschinen nach 
dem von HurtcHINSON und ScarroTT angegebenen Prinzip *)*). Bei Verwendung 
von Nickel machen allerdings die schlechte Reproduzierbarkeit der Grésse und 
Form der Empfangssignale sowie die betrachtliche Dampfung und Dispersion das 
Erreichen grosser Speicherkapazitaten schwierig. Gegen die Verwendung anderer 
Materialien spricht zunachst die in den meisten Fallen zu geringe Amplitude der 
empfangenen Signale. Eine andere Schwierigkeit besteht im Auftreten von Re- 
flexionssignalen bei ungeniigend gedampfter Befestigung der Leitungsenden. 

Zur Verbesserung der Betriebssicherheit und Erhéhung der Speicherkapazitat 
verwenden wir daher anstelle der iiblichen einzelnen Wandlerspule eine aus vielen 
Spulen bestehende Anordnung (Figur 1). 


SPEICHERDRAHT 


384 CU- WICKLUNG 
NR 56 261 10/228 
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(0) 10 20 30mm M1 FxD, wr. se 590. 76/15- ’ 


Figur 1 
Schnitt durch eine aus 12 Spulen bestehende elektroakustische Impulswandleranordnung. Bei Ver- 
wendung von Ni- oder Stahldrahten als akustische Leiter betragt die Zeitdifferenz zwischen auf- 
einanderfolgenden Spulen etwa 0,35 us. Die Vormagnetisierringe FXD 56.590.78/15C sind nur bei 
empfangsseitiger Anwendung erforderlich. i 


Dem Speicherdraht wird eine in der gewiinschten Richtung laufende elastische 
Welle aufgepragt, indem man die Spulen zeitlich nacheinander erregt. Die Lauf- 
zeitanpassung erfolgt durch Zusammenschalten mit Hilfsspulen und Kapazitaten 
gemass Figur 2. Mit je 12 Wandlerspulen erhielten wir sowohl sende- als auch — 
empfangsseitig je eine 6-10fache Erhdhung der Signalamplitude fiir die bevor- 
zugte Fortpflanzungsrichtung und nahezu vollstandige Ausléschung fiir die Ge- : 
genrichtung. Ausserdem stellten wir fest, dass sich die Form der Impulse leicht : 


1) Physikalisches Institut der ETH. ; 


2) G. W. Hurcuinson und G. G. Scarrort, Phil. Mag. [2] 42, 792 (1951): C. W. M 
, a Ge ; F : ; C. W. McCurTcHEn, 
Phil. Mag. [8] 2, 113 (1957). See “a 


3) D. Marper, Helv. phys. Acta 29, 459 (1956). 
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Figur 2 


Zusammenschaltung der Sende- bzw. Empfangsspulen S, ... Sy, (zum Beispiel je 100 Windungen, 

30 wH) mit den Hilfsspulen Sj... S{; und Kapazitaten (zum Beispiel je 1000 pF) zu einer Lauf- 

zeitkette (200 2/3,8 ws/0,6 ws). Ein am linken [rechten] Ende angeschlossener Sender pragt dem 
[bzw. Empfanger nimmt vom] Speicherdraht nur von links nach rechts laufende Impulse auf. 


durch Anderungen in der Art des Zusammenschaltens der Spulensysteme ver- 
bessern lasst. Der hohe Wirkungsgrad der neuen Wandleranordnung erlaubt nun 
auch die Verwendung von Stahldrahten mit dem Vorteil der geringeren Disper- 
sion und besseren Reproduzierbarkeit der Signalform. So erhielten wir nach einer 
Leitungslange von 18 m eines rostfreien Stahldrahtes (@ 0,15 mm) Signale von 
0,9 ws Basisbreite, was einer Speicherkapazitat von etwa 2000 Dualstellen ent- 
spricht. Bei Beniitzung dieser Leitung in einem dekadischen Addiersystem$) 


konnten wir 300 Dezimalstellen speichern. 


Fir die Eingangsverstarkung zwischen dem empfangsseitigen Wandler und 
der Addierschaltung hat sich ein Transistorenverstarker mit automatischer Pegel- 
regelung als niitzlich erwiesen. Ein ausfiihrlicher Bericht erscheint demndchst in 


den Helvetica Physica Acta. 


Herbst-Tagung der SPG 1957 


Die nachste Tagung der SPG (Schweizerische Physikalische Gesellschaft) 
wird am 22. September in Neuchatel stattfinden (Sekretariat: Basel, Klingel- 
bergstrasse 82). Die SPG wird als Sektion und im Rahmen der 137. Jahres- 
versammlung der Schweizerischen Naturforschenden Gesellschaft tagen (21. bis 


23. September, Neuchatel). 


K. P. MEYER 


Internationale Tagung fiir Elektronenmikroskopie 1958 in Berlin 


Die alle vier Jahre stattfindende Internationale Tagung der International 
Federation of Electron Microscope Societies wird im Jahre 1958 auf Einladung 
der Deutschen Gesellschaft fiir Elektronenmikroskopie in Berlin abgehalten. Das 
- Organisationskomitee leitet Prof. Ruska, Pradsident der International Federation 

of Electron Microscope Societies, Berlin-Dahlem, Faradayweg. Auskiinfte erteilt 
die Deutsche Gesellschaft fiir Elektronenmikroskopie, durch die Adresse Farb- 


werke Hoechst AG, Frankfurt am Main, zuhanden von Dr. KEHLER. 


7 : : P 


H. KEHLER 
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Buchbesprechungen — Book Reviews — Notices bibliographiques 


Journal of Fluid Mechanics. Band 1, Heft 1, Mai 1956. Herausgeber: 
Dr. G. K. BatcHEetor, Cavendish Laboratory, University of Cambridge, Cam- 
bridge, England. Druck und Verlag: Taylor & Francis Ltd., Red Lion Court, 
Fleet Street, London, E. C. 4. Preis: £1 ($3.00) pro Heft; Subskription Eppa 1s 0)010, 
pro Band mit sechs Heften. 

Das erste Heft dieser neuen Zeitschrift enthalt kein Vorwort, aus welchem die 
Absichten der Herausgeber hervorgingen. Die Zusammensetzung der Redaktion 
aus hervorragenden Theoretikern, wie G. K. BATCHELOR und M. J. LIGHTHILL, 
lasst immerhin vermuten, dass vor allem die physikalischen und mathematischen ~ 
Aspekte der Strémungslehre behandelt werden sollen, weniger aber die flug- 
technischen oder industriellen Anwendungen dieses Zweiges der Mechanik. Die — 
erste Nummer bestatigt diesen Eindruck. Unter anderem bringt sie Arbeiten tuber ~ 
die nichtlineare Theorie der Uberschallstrémung, eine interessante Untersuchung ~ 
der Sekundarstr6mung um Kugel und Zylinder, die experimentelle Bestimmung ~ 
der Relaxationszeiten fiir O,- und N,-Molekeln, eine Studie tiber die Grenzschicht — 
auf oszillierenden Platten ome 7 qindern sowie die Berechnung der Diffraktion 
von Stosswellen durch eine Kante. 

Wenn durch die Initiative der Herausgeber die bisher in verschiedenen Zeit- — 
schriften der physikalischen oder mathematischen Richtung zerstreuten Arbeiten — 
zusammengefasst und leichter zugdnglich gemacht werden k6énnen, so ist diese 
Neuerscheinung nur zu begriissen. P. de Hallev™ 


Jahrbuch der Wissenschaftlichen Gesellschaft fiir Luftfahrt e. V. 
(WGL) 1955. Herausgegeben von HERMANN BLENK (Verlag Friedrich a f 
& Sohn, Braunschweig: 1956). 360 S., 474 Abb.; DM 48.— 

Das vorliegende Jahrbuch der Wissenschaftlichen Gesellschaft fiir Luftfahrt 
e. V. gibt wiederum Zeugnis von der positiven, auf guter Tradition basierenden 
Tatigkeit dieser Gesellschaft. 1912 gegriindet, hat sie im Jahre 1936 beim Uber- | 
gang in die «Lilienthal-Gesellschaft fiir Luftfahrtforschung» eine Bilanz, beste-— 
hend aus 24 Jahrgangen der «Zeitschrift fiir Flugtechnik und Motorluftschiff-— 
fahrt» und aus 12 Jahrbiichern der WGL, aufzuweisen gehabt. Nachdem im 
Jahre 1945 die «Lilienthal-Gesellschaft» ihre Tatigkeit eingestellt hat, entstand 
in der deutschsprachigen, flugtechnischen Literatur eine Liicke, die aufzufiillen, 
unter anderen sich die 1952 wieder ins Leben gerufene WGL vorgenommen hat. — 
Als Mittel zur Erfiillung dieser Aufgabe dienen der Gesellschaft die « Zeitschrift 
fiir Flugwissenschaften (ZEW) » und die Jahrbticher der WGL. 

Der seit der Neugriindung 4, Jahrgang des Jahrbuches bringt vor allem Vor- 
trage der WGL-Tagung in Augsburg (12. bis 15. Oktober 1955) zum Abdruck, 
die, unter internationaler Beteiligung namhafter Fachleute, einen interessanten 
Querschnitt der flugtechnischen Probleme und Fragen der Gegenwart geben. 
Durch die Verdffentlichung sind Arbeiten aus folgenden Teilgebieten erfasst 
worden: Geschichte und Rechtsfragen, Tragfliigel-Aerodynamik und Gasdyna- 
mik, Propeller, Turbomaschinen und Strahltriebwerke, Betriebssicherheit, Steue- 
rung, Instrumente und Navigation, Messtechnik und Windkanale, Larmbeksaa 
fung, Flugzeugbau und Werkstoffe. 

Die wertvolle Ausgabe berechtigt, weitere Publikationen dieser regen Arbeits- 
gruppe mit aufmerksamem Interesse zu erwarten. Z. Plaskowshi 
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Prof, Dr. HELMUT Hasse 


Mathematik als Wissenschaft, 
Kunst und Macht 


32 Seiten, DINA5 Pappband mit Schutzumschlag DM 4.80 


Meine Ausfiihrungen richten sich keineswegs nur an Mathematiker, vielmehr an 
alle geistig und kulturell interessierten Menschen, mogen sie der mathematischen 
Wissenschaft nahe- oder fernstehen, sie hoch oder gering schatzen, wenn sie nur das 
Bediirfnis haben, ihren Blick einmal tiber ihr eigenes fachliches Tatigkeitsgebiet 
hinaus zu erheben, um sich fiir eine Stunde mit einer so bedeutenden Kulturerschei- 
nung zu beschaftigen, wie es die mathematische Wissenschaft ist. 

Nicht zuletzt wende ich mich auch an die akademische Jugend meines eigenen 
Faches; an diese Jugend, die im nur allzu voll ausgefiillten Rahmen des eigentlichen 
Fachstudiums erfahrungsgemass nur wenig Zeit findet, um sich tiber den tieferen 
Sinn der Mathematik und ihre Stellung im Ganzen unserer Kultur Gedanken zu 
machen. ; 

Eine Wissenschaft, deren iiber zweitausendjahrige Geschichte ein dauerndes, stetiges 
Aufsteigen in immer gréssere Hdéhen gewesen ist, darf stolz und mutig vor die Welt 
hintreten mit dem Bekenntnis: «Hier bin ich; seht voll Ehrfurcht auf die tiefen 
Erkenntnisse im Reiche des Denkens, die ich gewonnen habe, bewundert meine 
kristallklare, lebendige Schénheit und staunt, welche Fahigkeiten ich euch in die 
Hand gebe!» 
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Die erste Auflage wurde in wissenschaftlichen und technischen Zeitschriften sehr giinstig besprochen. 
Das Buch willin einem kurzen Uberblick die Kenntnisse vermitteln, die notwendig sind, um den 
Verbleib der Warmeenergie in technischen Betrieben tiberwachen und ausnutzen zu kénnen. Es gibt 


die Grundlagen der technischen Warmelehre in engstem Anschluss an die Praxis und unterscheidet — 


sich gerade durch diese unmittelbare Bezugnahme auf in der Technik auftretende Probleme von 
anderen Biichern, die diesen Grundlagen gewidmet sind. Der der 2. Auflage zugefiigte Anhang ent- 
halt teils Erganzungen fiir Leser, die an den inneren physikalischen Zusammenhangen besonderes 
Interesse haben, teils Hinweise auf Neuerungen, teils neue Beispiele. Der Text selbst wurde trotz 
der weitgehenden Verwendung des anastatischen Druckverfahrens an vielen Stellen wesentlich ver- 


andert, um neuen Erkenntnissen Rechnung zu tragen; auch neue Abbildungen wurden eingefiigt. 


In einer knappen, gut gegliederten Form und durch die tibersichtliche, das Wesentliche betonende 
Darstellungsweise ist das Buch den in der Praxis stehenden Ingenieuren ein schneller Ratgeber, 
insbesondere aber auch heranwachsenden Ingenieuren, z.B. den Studenten an den. technischen 


. Hochschulen und Universitaten, die ein technisches Fach ergriffen haben und sich-mit warmewirt- 


schaftlichen Fragen beschaftigen oder sich ganz allgemein der technischen und der experimentellen 
Physik widmen, beim Studium und zur Wiederholung von grossem Nutzen. - i 
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